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1 Introduction 



The standard big bang theory asserts that the Universe was born at some moment t — about 15 
billion years ago, in a state of infinitely large density and temperature. With the rapid expansion 
of the Universe the average energy of particles, given by the temperature, decreased rapidly, 
and the Universe became cold. This theory became especially popular after the discovery of 
the microwave background radiation. However, by the end of the 70's it was understood that 
this theory is hardly compatible with the present theory of elementary particles (primordial 
monopole problem, Polonyi fields problem, gravitino problem, domain wall problem) and it has 
many internal difficulties (flatness problem, horizon problem, homogeneity and isotropy problems, 
etc.). 

Fortunately, all these problems can be solved simultaneously in the context of a relatively 
simple scenario of the Universe evolution — the inflationary Universe scenario |l|j-f|. The main 
idea of this scenario is that the Universe at the very early stages of its evolution expanded 
quasi-exponentially (the stage of inflation) in a state with energy density dominated by the 
potential energy density V(<p) of some scalar field <fi. This rapid expansion made the Universe 
flat, homogeneous and isotropic, and made the density of monopoles, gravitinos and domain walls 
vanishingly small. Later, the potential energy density of the scalar field transformed into thermal 
energy, and still later, the Universe was correctly described by the standard hot Universe theory 
predicting the existence of the microwave background radiation. 

The first models of inflation were formulated in the context of the big bang theory. Their 
success in solving internal problems of this theory apparently removed the last doubts concerning 
the big bang cosmology. It remained almost unnoticed that during the last ten years inflationary 
theory changed considerably. It has broken an umbilical cord connecting it with the old big bang 
theory, and acquired an independent life of its own. For the practical purposes of describing the 
observable part of our Universe one may still speak about the big bang, just as one can still 
use Newtonian gravity theory to describe the Solar system with very high precision. However, 
if one tries to understand the beginning of the Universe, or its end, or its global structure, then 
some of the notions of the big bang theory become inadequate. For example, one of the main 
principles of the big bang theory is the homogeneity of the Universe. The assertion of homogeneity 
seemed to be so important that it was called "the cosmological principle" 0. Indeed, without 
using this principle one could not prove that the whole Universe appeared at a single moment of 
time, which was associated with the big bang. So far, inflation remains the only theory which 
explains why the observable part of the Universe is almost homogeneous. However, almost all 
versions of inflationary cosmology predict that on a much larger scale the Universe should be 
extremely inhomogeneous, with energy density varying from the Planck density to almost zero. 
Instead of one single big bang producing a single-bubble Universe, we are speaking now about 
inflationary bubbles producing new bubbles, producing new bubbles, ad infinitum [[U], |llj . Thus, 
recent development of inflationary theory considerably modified our cosmological paradigm ||. 
In order to understand better this modification, we should remember the main turning points in 
the evolution of the inflationary theory. 
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The history of inflationary cosmology goes back to 1965, to the papers by Erast Gliner [|IJ. He 
suggested that the Universe should begin its expansion in the vacuum-like state, and even calcu- 
lated the desirable amount of inflation ~ e 70 . However, at that time almost nobody took his ideas 
seriously, except for Andrei Sakharov, who made an attempt to calculate density perturbations 
produced in this scenario [|J. 



The first semi-realistic version of inflationary cosmology was suggested by Alexei Starobinsky 
Q. However, originally it was not quite clear what should be the initial state of the Universe 
in this scenario. Inflation in this model could not occur if the Universe was hot from the very 
beginning. To solve this problem, Zeldovich in 1981 suggested that the inflationary Starobinsky 
Universe was created "from nothing" [JT2 . This idea, which is very popular now [13|], at that time 



seemed too extravagant, and most cosmologists preferred to study inflation in a more traditional 
context of the hot Universe theory. 



One of the most important stages of the development of the inflationary cosmology was related 
to the old inflationary Universe scenario by Alan Guth Q . His scenario was based on the theory 
of the first order cosmological phase transitions, which was developed in the mid-70's by David 
Kirzhnits and myself |~4|]. This scenario was based on three fundamental propositions: 



1. The Universe initially expands in a state with a very high temperature, which leads to 
the symmetry restoration in the early Universe, <fi(T) = 0, where <ft is some scalar field driving 
inflation (the inflaton field). 

2. The effective potential V((j), T) of the scalar field has a deep local minimum at = 
even at a very low temperature T. As a result, the Universe remains in a supercooled vacuum 
state = (false vacuum) for a long time. The energy-momentum tensor of such a state rapidly 
becomes equal to T^ u = g^V^O), and the Universe expands exponentially (inflates) until the false 
vacuum decays. 



3. The decay of the false vacuum proceeds by forming bubbles containing the field 0o corre- 
sponding to the minimum of the effective potential V(4>). Reheating of the Universe occurs due 
to the bubble-wall collisions. 



The main idea of the old inflationary Universe scenario was very simple and attractive, and the 
role of the old inflationary scenario in the development of modern cosmology was very important. 
Unfortunately, as it was pointed out by Guth in [[|, this scenario had a major problem. If 
the rate of the bubble formation is bigger than the speed of the Universe expansion, then the 
phase transition occurs very rapidly and inflation does not take place. On the other hand, if the 
vacuum decay rate is small, then the Universe after the phase transition becomes unacceptably 
inhomogeneous. 

All attempts to suggest a successful inflationary Universe scenario failed until cosmologists 
managed to surmount a certain psychological barrier and renounce the aforementioned assump- 
tions, while retaining the main idea of ref. Q that the Universe has undergone inflation during the 
early stages of its evolution. The invention of the new inflationary Universe scenario || marked 



3 



the departure from the assumptions (2), (3). Later it was shown that the assumption (1) also 
does not hold in all realistic models known so far, for two main reasons. First of all, the time 
which is necessary for the field to roll down to the minimum of V((f), T) is typically too large, so 
that either inflation occurs before the field rolls to the minimum of V{<f>, T), or it does not occur 
at all. On the other hand, even if the field occasionally was near the minimum of V((j), T) from 
the very beginning, inflation typically starts very late, when thermal energy drops from M 4 = 1 
down to V(0,T). In all realistic models of inflation 1/(0, T) < 10~ 10 , hence inflation may start 
in a state with = not earlier than at t ~ 10 4 . During such a time a typical closed Universe 
would collapse before the conditions necessary for inflation could be realized f|. 

The assumption (1) was finally given up with the invention of the chaotic inflation scenario ||. 
The main idea of this scenario was to abandon the assumption that the Universe from the very 
beginning was hot, and that the initial state of the scalar field should correspond to a minimum 
of its effective potential. Instead of that, one should study various initial distributions of the 
scalar field 0, including those which describe the scalar field outside of its equilibrium state, and 
investigate in which case the inflationary regime may occur. 



In other words, the main idea of chaotic inflation was to remove all unnecessary restrictions 
on inflationary models inherited from the old theory of the hot big bang. In fact, the first step 
towards this liberation was already made when it was suggested to consider quantum creation 
of inflationary Universe from nothing [T2|, [Tj|. Chaotic inflation naturally incorporates this idea 
13], but it is more general: inflation in this scenario can also appear under a more traditional 



assumption of initial cosmological singularity. 



Thus, the main idea of chaotic inflation is very simple and general. One should study all 
possible initial conditions without insisting that the Universe was in a state of thermal equilib- 
rium, and that the field was in the minimum of its effective potential from the very beginning. 
However, this idea strongly deviated from the standard lore and was psychologically difficult to 
accept. Therefore, for many years almost all inflationary model builders continued the investi- 
gation of the new inflationary scenario and calculated high-temperature effective potentials in 
all theories available. It was argued that every inflationary model should satisfy the so-called 



'thermal constraints' [15|, that chaotic inflation requires unnatural initial conditions, etc. 



At present the situation is quite opposite. If anybody ever discusses the possibility that 
inflation is initiated by high-temperature effects, then typically the purpose of such a discussion 
is to show over and over again that this idea does not work (see e.g. |16|)) even though some 
exceptions from this rule might still be possible. On the other hand, there are many theories 
where one can have chaotic inflation (for a review see ||). 

A particularly simple realization of this scenario can be achieved in the theory of a massive 
scalar field with the effective potential '^(p 2 , or in the theory ^0 4 , or in any other theory with an 
effective potential which grows as n at large (whether or not there is a spontaneous symmetry 
breaking at small 0). Therefore a lot of work illustrating the basic principles of chaotic inflation 
was carried out in the context of these simple models. However, it would be absolutely incorrect 
to identify the chaotic inflation scenario with these models, just as it would be incorrect to identify 
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new inflation with the Coleman- Weinberg theory. The dividing line between the new inflation and 
the chaotic inflation was not in the choice of a specific class of potentials, but in the abandoning 
of the idea that the high-temperature phase transitions should be a necessary pre-requisite of 
inflation. 



Indeed, already in the first paper where the chaotic inflation was proposed ||, it was em- 
phasized that this scenario can be implemented not only in the theories ~ n , but in any model 
where the effective potential is sufficiently flat at some 0. In the second paper on chaotic inflation 



TTfl this scenario was implemented in a model with an effective potential of the same type as 



those used in the new inflationary scenario. It was explained that the standard scenario based on 
high-temperature phase transitions cannot be realized in this theory, whereas the chaotic inflation 
can occur there, either due to the rolling of the field from > 1, or due to the rolling down from 
the local maximum of the effective potential at = 0. In JTS| it was pointed out that chaotic 
inflation can be implemented in many theories including the theories with exponential potentials 
~ e°"^ with a < \Zl6n. This is precisely the same class of theories which was used to obtain the 
power law inflation . The class of models where chaotic inflation can be realized includes also 
the models based on the SU(5) grand unified theory ||, the R 2 inflation (modified Starobinsky 
model) POfl , supergravity-inspired models with polynomial and non-polynomial potentials pTf , 
'natural inflation' p2|, 'extended inflation' M, fl23|, 'hybrid inflation' ||24j| , etc. 



Thus, from the point of view of initial conditions, all inflationary models considered at present 
belong to the general class of chaotic inflation. We will describe the simplest version of chaotic 
inflation theory in Section 2, and in Section 3 we will discuss the problem of initial conditions for 
inflation. 



However, the way inflation begins is not the only possible way to classify inflationary models. 
Another classification describes various regimes which are possible during inflation: quasiexpo- 
nential inflation, power law inflation, etc. This classification is absolutely independent of the issue 
of initial conditions. Therefore it does not make any sense to compare, say, power law inflation 
and chaotic inflation, and to oppose them to each other. For example, in |1| it was pointed out 
that the chaotic inflation scenario, as distinct from the new inflationary Universe scenario, can 
be realized in the theories with the effective potential e a< ^ for a <C VTEn. Meanwhile, in [19] it 



was shown that this inflation is power law. Thus, the inflationary Universe scenario in the theory 
e a<p describes chaotic power law inflation. 



Finally, the third classification is related to the way inflation ends. There are three different 
possibilities: slow rollover, the first-order phase transition and the waterfall. The models of the 
first class describe slow rolling of the inflaton field 0, which gradually becomes faster and faster. 
A particular model of this type is chaotic inflation in the theories n . The models of the second 
class should contain at least two scalar fields, and a. They describe a strongly first-order phase 
transition with bubble production which is triggered by the slow rolling of the field 0. One of the 
popular models of this type is the extended inflation scenario 0, which is a combination of the 
Brans-Dicke theory and the old inflationary scenario. 

In the beginning it was assumed that the bubbles formed during the first-order phase transition 
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could be a useful ingredient of the theory of the large scale structure formation. However, later 
it was realized that one should make considerable modifications of the original models in order 
to avoid disastrous consequences of the bubble production. One should consider Brans-Dicke 
theories with a specifically adjusted effective potential and with a non-minimal kinetic term 
for the Brans-Dicke field, which does not seem pretty natural. According to the most recent 



modification |p5 ], the bubble formation happens only after the end of inflation. In this case, the 



end of inflation occurs as in the standard slow-rollover scenario. 

The waterfall regime appears in the models of two interacting scalar fields where inflation 



may end by a rapid rolling of the field a triggered by the slow rolling of the field [23|. This 
regime differs both from the slow-rollover and the first-order inflation. By changing parameters 
of this model one can continuously interpolate between these two regimes. Therefore we called 
this model "hybrid inflation". Hybrid models of such type may share the best features of the 
slow-rollover and the first-order models. We will describe the hybrid inflation model in Section 
4. 

One of the most important part of inflationary cosmology is the theory of reheating of the 
Universe after inflation. In fact, the word "reheating" is somewhat misleading since in most of 
the inflationary models which we study at present the assumption that the Universe was hot 
before inflation is unnecessary, at typically it is even harmful. Nevertheless, we are not going to 
change the old terminology at this point. It is much more important that the whole theory of 
reheating recently has been considerably changed. Previously we thought that reheating could be 
described by the standard methods of the theory of particle decay being applied to the inflaton 
scalar field driving inflation f26fl . Recently it was understood that in many theories the process 



of reheating is much more complicated, and it typically begins with a stage of an extremely rapid 
decay of the coherently oscillating scalar field due to a broad parametric resonance |2^j. Later 



this resonance becomes narrow, and finally the old theory of particle decay becomes applicable. 
However, this old theory should be applied not to the decay of the inflaton field, but to the decay 
of the particles produced at the first stages of explosive reheating, which we called "pre-heating" . 
The corresponding theory will be described in Section 5. 

One of the main purposes of inflationary theory was to solve the primordial monopole prob- 
lem. The solution was that inflation makes the distance between different monopoles exponen- 
tially large. However, recently we have learned that during inflation monopoles themselves could 
inflate. Moreover, inflation of monopoles continues without end even when inflation ends in the 
surrounding space. Therefore primordial monopoles (as well as other topological defects produced 
during inflation) can serve as seeds for the process of eternal self-reproduction of inflationary uni- 
28fl . This issue will be discussed in Section 6. 



verse 



Eternal inflation supported by monopoles is a very interesting and unexpected effect. However, 
eternal inflation is a very general property of inflationary models, which does not require any 
topological defects for its realization. In particular, this effect exists in the simplest versions of 



the chaotic inflation scenario [10], [ill : If the Universe contains at least one inflationary domain 
of a size of horizon (/i-region) with a sufficiently large and homogeneous scalar field 0, then 
this domain will permanently produce new /i-regions of a similar type. During this process the 
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total physical volume of the inflationary Universe (which is proportional to the total number of 
/i-regions) will grow indefinitely. The process of self-reproduction of inflationary domains occurs 
not only in the theories with the effective potentials growing at large <fi |IU], P] , but in some 
theories with the effective potentials used in old, new and extended inflation scenarios as well 
29| , |30| , EH) , j23| . However, in the models with the potentials growing at large the existence of 



this effect was most unexpected, and it leads to especially interesting cosmological consequences 



Self-reproduction of inflationary domains is responsible for the fundamental stationarity which 
is present in many inflationary models: properties of the parts of the Universe formed in the 
process of self-reproduction do not depend on the time when this process occurs. We call this 
property of the inflationary Universe local stationarity. In addition to it, there may exist a 
stationary distribution of probability P p to find a given field at a given time in a given physical 
volume. In this case we will be speaking of a global stationarity of the inflationary Universe. 



Inflationary models are divided into several different groups with respect to the self-reproduction 
and stationarity. In some of them inflation is possible but there is no self-reproduction. In others 
self- reproduction and the local stationarity are present, but there is no global stationarity. Fi- 
nally, there exist many models where the distribution P p is stationary ||32||-|p4||. We will consider 
this issue in Section 7. 



The possibility that our Universe is stationary (either locally or globally) shows that instead 
of being a fireball created 15 billion years ago, the Universe may be a huge expanding and self- 
reproducing fractal consisting of infinite number of exponentially large domains containing matter 
in all its possible phases. This changes our standard notions about what is natural and what is 
not. The new cosmological paradigm appears to be much more interesting and complicated than 
what we have expected. In the Section 8 we will briefly discuss some unusual nonperturbative 
effects which may imply that we should live in a center of a bubble, so that the density around 
us is smaller than the density at a very large distance from us (i.e. smaller than the critical 
density). This possibility is similar (though not quite identical) to the possibility that we live in 
an inflationary universe with Q < 1 |35| |. 

In all parts of this paper except for the Section 5 describing reheating we will use the system 
of units in which the Planck mass M p = 1. 



When I was finishing writing this paper I was informed about the death of academician 
Moisey Alexandrovich Markov. He was an extremely talented physicist who made a lot for the 
development of particle physics and theory of gravity, and also for breaking the artificial isolation 
of the Soviet physics from the international scientific community. I dedicate this article to his 
memory. 
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2 The simplest model of inflation 



We will begin our discussion with the simplest chaotic inflation model || based on the theory of 
a scalar field <fi minimally coupled to gravity, with the Lagrangian 

L = ^-R + \d,<pd^-V(<P) . (1) 

Here G = M~ 2 is the gravitational constant, R is the curvature scalar, and V((f>) is the effective 
potential of the scalar field. If the classical field <fi is sufficiently homogeneous in some domain of 
the Universe, then its behavior inside this domain is governed by the equations 

<j> + 3H<j) = -dV/d(f> , (2) 

Here H = a/a, a(t) is the scale factor of the Universe, k = +1, —1, or for a closed, open or flat 
Universe, respectively. 

The simplest version of the theory (jl]) is the theory of a massive noninteracting scalar field 
with the effective potential V(<f>) = \m 2 (j> 2 , where m is the mass of the scalar field <fi, m <C 1. If 
the field initially is sufficiently large, <fi ^> 1, then one can show that the functions <p(t) and 
a(t) rapidly approach the asymptotic regime 

<P(t) = O - ' fl W = a ° • ex P ( 2n " 2 W)) • ( 4 ) 

Note that in this regime the second derivative of the scalar field in eq. (0) and the term k/a 2 in 
(^) can be neglected. The last statement means that the Universe becomes locally flat. 

According to (^), during the time r ~ <p/m (which is much greater than H for <fi ^> 1) the 
relative change of the field remains small, the effective potential V{4>) changes very slowly 
and the Universe expands quasi-exponentially: a(t + At) ~ a(t) exp(HAt) for At < r. Here 
H = 2^lm(t). 

The regime of the slow rolling of the field and the quasi-exponential expansion (inflation) 
of the Universe ends at <fi < <p e . In the theory under consideration, <p e ~ 0.2. At <fi < <j) e the 
field <p oscillates rapidly, and if this field interacts with other matter fields (which are not written 
explicitly in eq. (P), its potential energy V(<f>) ~ HL^s. ^ ^_ j s transformed into heat. The 
reheating temperature Tr may be of the order m 1 ^ 2 or somewhat smaller, depending on the 
strength of the interaction of the field <fi with other fields. It is important that Tr does not 
depend on the initial value 0o of the field 0. The only parameter which depends on 0o is the 
scale factor a(t), which grows e 27r<?i ° times during inflation. 

All results obtained above can be easily generalized for the theories with more complicated 
effective potentials. For example, during inflation in the theories with V((j)) = - <p n one has 

^ (t ) = ^-izltg ( W 4, (5) 
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<j>(t) = O exp -J — t) for n = 4 . (6) 



Inflation ends at e ~ 4 JL- . The scale factor is given by 

a(t) = a o exp(^(0 2 ,- 2 (t))) . (7) 



Note, that in all realistic models of elementary particles spontaneous symmetry breaking 
occurs on a scale which is many orders of magnitude smaller than 1 (i.e. smaller than M p ). 
Therefore all results which we obtained remain valid in all theories which have potentials V{4>) ~ 
- 4> n at > 1, independently of the issue of spontaneous symmetry breaking which may occur in 
such theories at small 0. 

As was pointed out in |19|, [0| , chaotic inflation occurs as well in the theories with exponential 



effective potentials V(0) = Vq e a ^ for sufficiently small a. In such theories the Einstein equations 
and equations for the scalar field have an exact solution: 

0(t) = 0o-- , a(t) = a t p , (8) 

a to 

where p = ^f. Note that here we are dealing with the power law expansion of the Universe. It 
can be called inflation if p ^> 1, which implies that a <C \A67r ~ 7. 



3 The way inflation begins 



It is most important to investigate initial conditions which are required for inflation. Let us 
consider first a closed Universe of initial size I ~ 1 (in Planck units), which emerges from the 
space-time foam, or from singularity, or from 'nothing' in a state with the Planck density p ~ 1. 
Only starting from this moment, i.e. at p < 1, can we describe this domain as a classical Universe. 
Thus, at this initial moment the sum of the kinetic energy density, gradient energy density, and 
the potential energy density is of the order unity: |0 2 + |(<9j0) 2 + V(0) ~ 1. 

We wish to emphasize, that there are no a priori constraints on the initial value of the scalar 
field in this domain, except for the constraint |0 2 + |(<9i0) 2 + V"(0) ~ 1. Let us consider for 
a moment a theory with V((p) = const. This theory is invariant under the shift — >• + a. 
Therefore, in such a theory all initial values of the homogeneous component of the scalar field 
are equally probable. Note, that this expectation would be incorrect if the scalar field should 
vanish at the boundaries of the original domain. Then the constraint |(<9j0) 2 < 1 would tell us 
that the scalar field cannot be greater than 1 inside a domain of initial size 1. However, if the 
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original domain is a closed Universe, then it has no boundaries. (We will discuss a more general 
case shortly.) 



The only constraint on the average amplitude of the field appears if the effective potential is 
not constant, but grows and becomes greater than the Planck density at > 4> p , where V(<p p ) = 1. 
This constraint implies that <p < <p p , but it does not give any reason to expect that (f) <C <p p . This 
suggests that the typical initial value O of the field in such a theory is </> ~ <p p . 

Thus, we expect that typical initial conditions correspond to |0 2 ~ ^(di<f)) 2 ~ V(<j>) = 0(1). 
Note that if by any chance |0 2 + |(<9j0) 2 < V(4>) in the domain under consideration, then inflation 
begins, and within the Planck time the terms |0 2 and |(<9j0) 2 become much smaller than V((f>), 
which ensures continuation of inflation. It seems therefore that chaotic inflation occurs under 
rather natural initial conditions, if it can begin at V(<p) ~ 1 || [TB[. 

The assumption that inflation may begin at a very large <fi has important implications. For 
example, in the theory @) one has 0o ~ 0p ~ m -1 / 2 . Let us consider for definiteness a closed 
Universe of a typical initial size 0(1). Then, according to (^), the total size of the Universe after 
inflation becomes equal to 

/ ~ exp (2tt0 2 ) ~ exp (-^) . (9) 

For m ~ 10" 6 (which is necessary to produce density perturbations — ~ 10 5 , see below) 

I ~ exp (27rl0 12 ) > 10 1C|12 cm. Thus, according to this estimate, the smallest possible domain of 
the Universe of initial size 0(M~ l ) ~ 10~ 33 cm after inflation becomes much larger than the size 
of the observable part of the Universe ~ 10 28 cm. This is the reason why our part of the Universe 
looks flat, homogeneous and isotropic. 

In what follows we will return many times to our conclusion that the most probable initial 
value of the scalar field corresponds to <fi ~ <p p . There were many objections to it. Even though 
all these objections were answered many years ago []8], [T^| , we need to discuss here one of these 



objections again. This is important for a proper understanding of the new picture of the evolution 
of the Universe in inflationary cosmology. 

Assume that the Universe is not closed but infinite, or at least extremely large from the very 
beginning. (This objection does not apply to the closed Universe scenario discussed above.) In 



this case one could argue that our expectation that 0o ~ 0p 3> 1 is not very natural |36j . Indeed, 
the conditions (<9j0) 2 < 1 and O ~ 4>p ^ 1 imply that the field (p should be of the same order of 
magnitude ~ <p p 3> 1 on a length scale at least as large as P , which is much larger than the 
scale of horizon / ~ 1 at the Planck time. But this is highly improbable, since initially (i.e., at 
the Planck time ) there should be no correlation between values of the field <fi in different regions 
of the Universe separated from one another by distances greater than 1. The existence of such 
correlation would violate causality. As it is written in [3B], the scalar field </> must be smooth on 



a scale much greater than the scale of the horizon, which does not sound very chaotic. 

The answer to this objection is very simple || |T8 |. We have absolutely no reason to expect 
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that the overall energy density p simultaneously becomes smaller than the Planck energy density 
in all causally disconnected regions of an infinite Universe, since that would imply the existence 
of an acausal correlation between values of p in different domains of Planckian size l p ~ 1. Thus, 
each such domain at the Planck time after its creation looks like an isolated island of classical 
space-time, which emerges from the space-time foam independently of other such islands. During 
inflation, each of these islands independently acquires a size many orders of magnitude larger 
than the size of the observable part of the Universe. A typical initial size of a domain of classical 
space-time with p < 1 is of the order of the Planck length. Outside each of these domains 
the condition p < 1 no longer holds, and there is no correlation between values of the field <j) 
in different disconnected regions of classical space-time of size 1. But such correlation is not 
necessary at all for the realization of the inflationary Universe scenario: according to the 'no hair' 
theorem for de Sitter space, a sufficient condition for the existence of an inflationary region of 
the Universe is that inflation takes place inside a region whose size is of order if -1 . In our case 
this condition is satisfied. 



We wish to emphasize once again that the confusion discussed above, involving the correlation 
between values of the field <fi in different causally disconnected regions of the Universe, is rooted in 
the familiar notion of a very large Universe that is instantaneously created from a singular state 
with p = oo, and instantaneously passes through a state with the Planck density p — 1. The lack 
of justification for such a notion is the very essence of the horizon problem. Now, having disposed 
of the horizon problem with the aid of the inflationary Universe scenario, we can possibly manage 
to familiarize ourselves with a different picture of the Universe. In this picture the simultaneous 
creation of the whole Universe is possible only if its initial size is of the order 1, in which case no 
long-range correlations appear. Initial conditions should be formulated at the Planck time and 
on the Planck scale. Within each Planck-size island of the classical space-time, the initial spatial 
distribution of the scalar field cannot be very irregular due to the constraint (<9j0) 2 < 1. But 
this does not impose any constraints on the average values of the scalar field <fi in each of such 
domains. One should examine all possible values of the field and check whether they could lead 
to inflation. 



Now let us analyze the probability for inflation to occur in the theories with small V(<p) . 
There are two for such an event to happen. First of all, the inflationary universe can appear from 
nothing in a state with corresponding to the maximum of V(<p). The probability of such an 
event can be estimated by the square of the tunneling wave function of the Universe [13] : 



p ~^(-wm)- (10) 

This probability may be quite high if inflation is possible for V((p) ~ 1, as in the chaotic inflation 
scenario. However, in the models where inflation occurs near the maximum of V(<fr) considerations 
related to the amplitude of density perturbations require V((f)) < 10~ 10 , which gives a rather 
disappointing result for the probability of inflation, P ~ exp(— 10 10 ). 

Another possibility is that the Universe was large and hot from the very beginning, and it 
became inflationary only when its thermal energy density ~ T 4 dropped below V(<p). Inflation 
begins at this stage only if the Universe at this moment was sufficiently homogeneous on a scale 
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greater than Ax ~ H ~ l/yV(0). Suppose for simplicity that the Universe from the very 
beginning was dominated by ultrarelativistic matter. Then its scale factor expanded as p -4 , 
where p is the energy density at the pre-inflationary stage. Therefore at the Planck time the 



size of the part of the Universe which later evolved into inflationary domain was not l/yV(<p), 
but somewhat smaller: Ax ~ V r_1 / 4 (0). This whole scenario can work only if at the Planck 
time the domain of this size was sufficiently homogeneous, ^ < 1. However, at the Planck time 

this domain consisted of V~ 3 ^((j)) domains of a Planck size, and energy density in each of them 
was absolutely uncorrelated with the energy density in other domains. Therefore a priori one 
could expect changes of density 5p ~ p when going from one causally disconnected parts of the 
Universe of a size Mp 1 = 1 to another. Simple combinatorial analysis suggests that the probability 
of formation of a reasonably homogeneous part of the Universe of a size Ax ~ V A_1/ ' 4 (0) at the 
Planck time is suppressed by the exponential factor 

where C = 0(1). To get a numerical estimate, one can take V(0) ~ 10 -10 . This gives P < 10^ 1q7 . 



Note, that the arguments given above || [TS[] suggest that initial conditions for inflation are 
quite natural only if inflation begins as close as possible to the Planck density. These arguments 
do not give any support to the models where inflation is possible only at densities much smaller 
than 1. Inflationary models of that type require fine-tuned initial conditions, and they cannot 
solve the flatness problem, unless we apply some rather sophisticated arguments based on the 



theory of a self-reproducing inflationary universe . 



4 The way inflation ends. Hybrid Inflation 

As we already mentioned in the introduction, there exist three different way for inflation to end. 
The first one is by slow rolling, which gradually becomes faster and faster. This possibility is 
illustrated by the simplest model considered in the second section. The second possibility is the 
first order phase transition 0. However, after many modifications of this scenario, it gradually 
became rather complicated. It requires specific adjustments of the effective potential and of the 
kinetic terms for the Brans-Dicke field |2|| . The third possibility is the rapid rolling of one of the 
scalar field triggered by the slow rolling of another (waterfall) ||24|| . We will discuss this possibility 
here. 

The simplest model where the waterfall regime can appear is the theory of two scalar fields 
with the effective potential 

V(a, 0) = ^(M 2 - A<x 2 ) 2 + ^ + |>V . (12) 
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Theories of this type were considered in |38[ -|40[. The main difference between the models of 



refs. [pq| -|[40|| and our model is a specific choice of parameters, which allows the existence of the 
waterfall regime mentioned above. There is also another important difference: we will assume 
that the field a in this model is the Higgs field, which remains as a physical degree of freedom 
after the Higgs effect in an underlying gauge theory with spontaneous symmetry breaking. This 
field acquires only positive values, which removes the possibility of domain wall formation in this 
theory. Usually it is rather dangerous to take the inflaton field interacting with gauge fields, since 
its effective coupling constant A may acquire large radiative corrections ~ e 4 , where e is the gauge 
coupling constant. In our case this problem does not appear since density perturbations in our 
model remain small at rather large A, see below. 

The effective mass squared of the field a is equal to —M 2 + g 2 (p 2 . Therefore for > C = M/g 
the only minimum of the effective potential V(a, 0) is at a = 0. The curvature of the effective 
potential in the cr-direction is much greater than in the 0-direction. Thus we expect that at the 
first stages of expansion of the Universe the field a rolled down to o = 0, whereas the field 
could remain large for a much longer time. For this reason we will consider the stage of inflation 
at large 0, with o = 0. 

At the moment when the inflaton field becomes smaller than C = M/g, the phase transition 
with the symmetry breaking occurs. If m 2 cf) 2 c = m 2 M 2 /g 2 -C M 4 /A, the Hubble constant at the 
time of the phase transition is given by 

9 2vrM 4 , . 

H = — ■ (13) 

Thus we will assume that M 2 We will assume also that m 2 <C H 2 , which gives 

/ q\ 

M 2 »mJ- . (14) 
V 2n 

One can easily verify, that, under this condition, the Universe at > C undergoes a stage of 
inflation. In fact, inflation in this model occurs even if m 2 is somewhat greater than H 2 . Note 
that inflation at its last stages is driven not by the energy density of the inflaton field but by 
the vacuum energy density V(0, 0) = as in the new inflationary Universe scenario. This was 
the reason why we called this model 'hybrid inflation' in . 



Let us study the behavior of the fields and a after the time At = H~ l = \J 2 ^ti 2 f rom the 
moment t c when the field becomes equal to C . The equation of motion of the field during 
inflation is 3H<fi = m 2 (f). Therefore during the time interval At = H^ 1 the field decreases from 
0c by A0 = YjjT = 2wg l M s • ^ ne absolute value of the negative effective mass squared — M 2 + g 2 (f) 2 
of the field a at that time becomes equal to 

MV) = ^ ■ (15) 

The value of M 2 (0) is much greater than H 2 for M 3 cAm. In this case the field a within the 
time At ~ H^ 1 rolls down to its minimum at cr(0) = M(0)/a/A, rapidly oscillates near it and 
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loses its energy due to the expansion of the Universe. However, the field cannot simply relax near 
this minimum, since the effective potential V((f>, u) at cr(0) has a nonvanishing partial derivative 

^ = ^ + ^' (16) 

One can easily check that the motion in this direction becomes very fast and the field <j) rolls to 
the minimum of its effective potential within the time much smaller than H' 1 if M 3 <C V\g m. 
Thus, under the specified conditions inflation ends in this theory almost instantaneously, as soon 
as the field reaches its critical value 4> c = M/g. 

The amplitude of adiabatic density perturbations produced in this theory can be estimated 
by standard methods II] and is given by 

(17) 



p 5 5m 2 (f) 

In the case m 2 <C H 2 the scalar field <fi does not change substantially during the last 60 e-foldings 
(i.e. during the interval At ~ QOH^ 1 ). In this case the amplitude of density perturbations 
practically does not depend on scale, and is given by 



5p 2y/6^M 5 
p 5AvAm 2 

The definition of ^ used in || corresponds to COBE data for ~ 5 • 10~ 5 . Dividing it by 
(14) with an account taken of fll8|) gives M 3 <C 5 • 10~ 5 \g~ 1 m. This means that the 'waterfall 
conditions' M 3 <Am and M 3 <C \f\g m automatically follow from the conditions m 2 -C H 2 
and ^ ~ 5 • 10~ 5 , unless the coupling constants A and g are extremely small. Therefore the 
waterfall regime is realized in this model for a wide variety of values of parameters m, M, A and 
g which lead to density perturbations ~ 5 • 10~ 5 . 

To give a particular example, let us take g 2 ~ A ~ 10 _1 , m ~ 10 2 GeV (electroweak scale). 

In this case all conditions mentioned above are satisfied and ^ ~ 5 ■ 10^ 5 for M ~ 1.3 ■ 10 11 GeV. 

p 

In particular, we have verified, by solving equations of motion for the fields <\> and a numerically, 
that inflation in this model ends up within the time At H^ 1 after the field <fi reaches its critical 
value 4> c = M/g. The value of the Hubble parameter at the end of inflation is given by H ~ 7- 10 3 
GeV. The smallness of the Hubble constant at the end of inflation makes it possible, in particular, 
to have a consistent scenario for Hybrid in inflationary cosmology even if the axion mass is much 
smaller than 10~ 5 eV pl |. This model has some other distinctive features. For example, the 



spectrum of perturbations generated in this model may look as a power-law spectrum rapidly 



decreasing at large wavelength I pT . 



Indeed, at the last stages of inflation (for ^ 3> the field behaves as 

/ m 2 (t — t r )\ , , 

<P = <p c ■ exp ( ^-^) , (19) 
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whereas the scale factor of the Universe grows exponentially, a ~ e . This leads to the following 
relation between the wavelength of perturbations I and the value of the scalar field <fi at the 

/ , \m 2 /3H 2 

moment when these perturbations were generated: <\> ~ <p c ( — J . In this case 



5p 2V6vr^M 



r5 



p 5Aa/A 



(20) 



which corresponds to the spectrum index n = 1 + |^ = 1 + -^z- Note that this spectrum index 
is greater than 1, which is a very unusual feature. For the values of m, M, A and g considered 
above, the deviation of n from 1 is vanishingly small (which is also very unusual). However, let 
us take, for example, A = g — 1, M — 10 15 GeV (grand unification scale), and m — 5 • 10 10 GeV. 
In this case the amplitude of perturbations at the end of inflation (<f) = <ft c ) is equal to 4 ■ 1CT 4 , 
n ~ 1.1, and the amplitude of the density perturbations drops to the desirable level ~ 5 ■ 1CT 5 
on the galaxy scale (l g ~ l c ■ e 50 ). One may easily obtain models with even much larger n, but 



this may be undesirable, since it may lead to formation of many small primordial black holes |p7 |. 



Note, that the decrease of ^ at large / is not unlimited. At > ^ the spectrum begins 
growing again. Thus, the spectrum has a minimum on a certain scale, corresponding to the 
minimum of expression flT7D. This complicated shape of the spectrum appears in a very natural 
way, without any need to design artificially bent potentials. 

As we have seen, coupling constants in our model can be reasonably large, and the range of 
possible values of masses m and M is extremely wide. Thus, our model is very versatile. One 
should make sure, however, that the small effective mass of the scalar field does not acquire 
large radiative corrections near = <f> c . Hopefully this can be done in supersymmetric theories 
with flat directions of the effective potential. 



5 Reheating after inflation 



The theory of reheating of the Universe after inflation is the most important application of 
the quantum theory of particle creation, since almost all matter constituting the Universe at 
the subsequent radiation-dominated stage was created during this process |§. At the stage of 
inflation all energy was concentrated in a classical slowly moving inflaton field (p. Soon after the 
end of inflation this field began to oscillate near the minimum of its effective potential. Gradually 
it produced many elementary particles, they interacted with each other and came to a state of 
thermal equilibrium with some temperature T r , which was called the reheating temperature. 



An elementary theory of reheating was first developed in [26] for the new inflationary scenario. 
Independently a theory of reheating in the R 2 inflation was constructed in ||4"2"f . Various aspects 

H H 



of this theory were further elaborated by many authors, see e.g. 



Still, a general scenario 
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of reheating was absent. In particular, reheating in the chaotic inflation theory remained almost 
unexplored. The present section contains results obtained recently in our work with Kofman 
and Starobinsky P7| . A more detailed presentation of our results is to be published in JI3 



We have found that the process of reheating typically consists of three different stages. At 
the first stage, which cannot be described by the elementary theory of reheating, the classical 
coherently oscillating inflaton field decays into massive bosons (in particular, into 0-particles) 
due to parametric resonance. In many models the resonance is very broad, and the process 
occurs extremely rapidly (explosively). Because of the Pauli exclusion principle, there is no 
explosive creation of fermions. To distinguish this stage from the stage of particle decay and 
thermalization, we will call it pre-heating. Bosons produced at that stage are far away from 
thermal equilibrium and typically have enormously large occupation numbers. The second stage 
is the decay of previously produced particles. This stage typically can be described by methods 
developed in p6f . However, these methods should be applied not to the decay of the original 
homogeneous inflaton field, but to the decay of particles and fields produced at the stage of 
explosive reheating. This considerably changes many features of the process, including the final 
value of the reheating temperature. The third stage is the stage of thermalization, which can be 
described by standard methods, see e.g. [§, p6| ; we will not consider it here. Sometimes this stage 
may occur simultaneously with the second one. In our investigation we have used the formalism of 
the time-dependent Bogoliubov transformations to find the density of created particles, n^(t). A 
detailed description of this theory will be given in [43j; here we will outline our main conclusions 
using a simple semiclassical approach. 



sec 



We should note that recently two other groups of authors have also studied related problems, 
jBL pTTfl. The corresponding papers are very useful and contain many interesting results, but 



some additional work is to be done in order to apply these results to the theory of reheating 



after inflation. For example, ref. ||6| was devoted to study of the stage of a narrow resonance 



without an account taken of backreaction of created particles. As we already mentioned, the 
first stages of reheating typically occur in the regime of a broad resonance. On the other hand, 
as it was shown above, the stage of a broad resonance ends because of the backreaction of 
created particles, and at the stage of a narrow resonance the effects of backreaction are also very 
important. Backreaction of created particles was investigated in a very detailed way in ref. j|7| . 



However, their investigation was performed neglecting expansion of the Universe. In the theory 



,2 



^-(j> + with m > the authors studied the strong coupling regime, A ^> 1, whereas in the 
chaotic inflation scenario which we are investigating now A ~ 1CT 13 . Therefore in what follows 
we will concentrate on the results obtained in [pTl] . 



We will consider a simple chaotic inflation scenario describing the classical inflaton scalar field 
4> with the effective potential V(<p) = ±^m^(p 2 + j4> A . Minus sign corresponds to spontaneous 
symmetry breaking <p — > <p + o with generation of a classical scalar field a — 7^=. The field (f) 

after inflation may decay into bosons \ an d fermions ip due to the interaction terms — \g 2 4> 2 X 2 
and —hipipcf). Here A, g and h are small coupling constants. In case of spontaneous symmetry 
breaking, the term — \g 2 <fi 2 x 2 gives rise to the term —g 2 a<f)x 2 - We will assume for simplicity 
that the bare masses of the fields x an d ip are very small, so that one can write rn x ((f)) = gcf), 
m^(4>) = \h(j)\. 
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Let us briefly recall the elementary theory of reheating ||. At > M p , we have a stage of 
inflation. This stage is supported by the friction-like term 3if0 in the equation of motion for 
the scalar field. Here H = a/a is the Hubble parameter, a(t) is the scale factor of the Universe. 
However, with a decrease of the field this term becomes less and less important, and inflation 
ends at < M p /2. After that the field begins oscillating near the minimum of V{4>). The 
amplitude of the oscillations gradually decreases because of expansion of the Universe, and also 
because of the energy transfer to particles created by the oscillating field. Elementary theory 
of reheating is based on the assumption that the classical oscillating scalar field 0(t) can be 
represented as a collection of scalar particles at rest. Then the rate of decrease of the energy of 
oscillations coincides with the decay rate of 0-particles. The rates of the processes — > XX an d 
— > ipip (for ^> 2m x , 2m^) are given by 

r(0 - xx) = , r(0 - w>) = ^ . (21) 

oirm^ 87T 

Reheating completes when the rate of expansion of the Universe given by the Hubble constant 
H = J^m% ~ t~ l becomes smaller than the total decay rate V = T(0 — > xx) + T(0 — > ip4>)- The 



reheating temperature can be estimated by T r ~ 0.1 JTM. 



It is interesting to note that in accordance with the elementary theory of reheating the ampli- 
tude squared of the oscillating scalar field decays exponentially, as e~ Tt . Phenomenologically, this 
can be described by adding the term T0 to the equation of motion of the scalar field. Unfortu- 
nately, many authors took this prescription too seriously and investigated the possibility that the 
term T0, just like the term 3H<p, can support inflation. We should emphasize f2?H , that adding 



the term T0 to the equation of motion is justified only at the stage of oscillations (i.e. after the 
end of inflation), and only for the description of the amplitude of oscillations of the scalar field, 
rather than for the description of the scalar field itself. Moreover, even at the stage of oscillations 
this description becomes incorrect as soon as the resonance effects become important. 

As we already mentioned, elementary theory of reheating can provide a qualitatively correct 
description of particle decay at the last stages of reheating. Moreover, this theory is always 
applicable if the inflaton field can decay into fermions only, with a small coupling constant 
h 2 <C m^/Mp. However, typically this theory is inapplicable to the description of the first stages 
of reheating, which makes the whole process quite different. In what follows we will develop the 
theory of the first stages of reheating. We will begin with the theory of a massive scalar field 
decaying into particles x> then we consider the theory ^0 4 , and finally we will discuss reheating 
in the theories with spontaneous symmetry breaking. 

We begin with the investigation of the simplest inflationary model with the effective potential 

-^0 2 . Suppose that this field only interacts with a light scalar field x ( m x ^ m <t>) ^ ue t° the term 
— \g 2 <t> 2 X 2 - The equation for quantum fluctuations of the field x with the physical momentum 
k/a(t) has the following form: 

Xk + SH Xk + ( J^y + 9 2 ® 2 sin 2 (?V)) Xk = , (22) 
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where k = yk 2 , and $ stands for the amplitude of oscillations of the field 0. As we shall 
see, the main contribution to x-P ar ticl e production is given by excitations of the field x with 
k/a ^> rritj,, which is much greater than H at the stage of oscillations. Therefore, in the first 
approximation we may neglect the expansion of the Universe, taking a(t) as a constant and 
omitting the term 3Hxk m fl22[)- Then the equation (^) describes an oscillator with a variable 
frequency ^f(t) = k 2 a~ 2 + g 2 $ 2 sin 2 (m^t). Particle production occurs due to a nonadiabatic 
change of this frequency. Equation (^) can be reduced to the well-known Mathieu equation: 

Xk + (A(k)-2qcos2z) X k = , (23) 

where A(k) = -^r^ + 2q, q = z = m^t, prime denotes differentiation with respect to 

z. An important property of solutions of the equation ( p3| ) is the existence of an exponential 
instability Xk oc exp(/j,^ z) within the set of resonance bands of frequencies Ak^ labeled by 
an integer index n. This instability corresponds to exponential growth of occupation numbers 
of quantum fluctuations n^(t) oc exp (2//j^ m^i) that may be interpreted as particle production. 
As one can show, near the line A = 2q there are regions in the first, the second and the higher 
instability bands where the unstable modes grow extremely rapidly, with //& ~ 0.2. We will show 
analytically in |E| that for q ^> 1 typically [ik ~ ~ 0.175 in the instability bands along the 

line A = 2q, but its maximal value is ln ( 1 + v/2 ) ~ 0.28. Creation of particles in the regime of a 
broad resonance (q > 1) with 2irfj 1 k = 0(1) is very different from that in the usually considered 
case of a narrow resonance (q <C 1), where 2iTfik < 1 In particular, it proceeds during a tiny 
part of each oscillation of the field when 1 — cos z ~ q^ 1 and the induced effective mass of 
the field x (which is determined by the condition m 2 = g 2 $ 2 /2) is less than m^. As a result, 
the number of particles grows exponentially within just a few oscillations of the field 0. This 
leads to an extremely rapid (explosive) decay of the classical scalar field 0. This regime occurs 
only if q > 7r _1 , i.e. for g$ > m^, so that gM p is the necessary condition for it. One 

can show that a typical energy E of a particle produced at this stage is determined by equation 
A — 2q ~ y/q, and is given by E ~ yj gm^Mp f43fl . 

Creation of x-particles leads to the two main effects: transfer of the energy from the homo- 
geneous field <p(t) to these particles and generation of the contribution to the effective mass of 
the field: wl e // = m ^> + 9 2 (x 2 )ren- The last term in the latter expression quickly becomes 
larger than mi. One should take both these effects into account when calculating backreac- 
tion of created particles on the process. As a result, the stage of the broad resonance creation 
ends up within the short time t ~ mT 1 \n(m^/g 5 M p ), when $ 2 ~ (x 2 ) and q = — becomes 

smaller than 1. At this time the energy density of produced particles ~ E 2 (x 2 ) ~ gm^M v ^ 2 
is of the same order as the original energy density ~ m^M 2 of the scalar field at the end of 
inflation. This gives the amplitude of oscillations at the end of the stage of the broad resonance 
particle creation: $ 2 ~ (x 2 ) ~ g~ x m^,M p <^ M 2 . Since E ^> m^, the effective equation of 
state of the whole system becomes p ~ e/3. Thus, explosive creation practically eliminates a 
prolonged intermediate matter-dominated stage after the end of inflation which was thought to 
be characteristic to many inflationary models. However, this does not mean that the process of 
reheating has been completed. Instead of x-particles in the thermal equilibrium with a typical 
energy E ~ T ~ (mMp) 1 ^ 2 , one has particles with a much smaller energy ~ (gm^Mp) 1 ^ 2 , but 
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with extremely large mean occupation numbers ~ g ^> 1. 

After that the Universe expands as a(t) oc y/i, and the scalar field continues its decay in 
the regime of the narrow resonance creation q « <C 1. As a result, decreases rather 

slowly, oc t~ 3 / 4 . This regime is very important because it makes the energy of the field 
much smaller than that of the x _ P ar ticles. One can show that the decay finally stops when 
the amplitude of oscillations $ becomes smaller than g~ x m$ [43]. This happens at the moment 



t ~ m^ 1 (gM p /m ( f,) 1 ^ 3 (in the case m < g 7 M p decay ends somewhat later, in the perturbative 
regime). The physical reason why the decay stops is rather general: decay of the particles in 
our model occurs due to its interaction with another 0-particle (interaction term is quadratic 
in and in x). When the field (or the number of 0-particles) becomes small, this process is 
inefficient. The scalar field can decay completely only if a single scalar 0-particle can decay into 
other particles, due to the processes — > XX or ^ ijjijj, see eq. (|2"T]). If there is no spontaneous 
symmetry breaking and no interactions with fermions in our model, such processes are impossible. 

At later stages the energy of oscillations of the inflaton field decreases as a~ 3 (t), i.e. more 
slowly than the decrease of energy of hot ultrarelativistic matter oc a~ 4 (t). Therefore, the relative 
contribution of the field 0(t) to the total energy density of the Universe rapidly grows. This gives 
rise to an unexpected possibility that the inflaton field by itself, or other scalar fields can be cold 
dark matter candidates, even if they strongly interact with each other. However, this possibility 
requires a certain degree of fine tuning; a more immediate application of our result is that it 
allows one to rule out a wide class of inflationary models which do not contain interaction terms 
of the type of g 2 a<px 2 or hetnpip. 

So far we have not considered the term ^0 4 in the effective potential. Meanwhile this term 
leads to production of 0-particles, which in some cases appears to be the leading effect. Let us 

study the 0-particle production in the theory with V(0) = ^r0 2 + |0 4 with m 2 «C AM 2 . In this 
case the effective potential of the field soon after the end of inflation at ~ M p is dominated 
by the term ^0 4 . Oscillations of the field in this theory are not sinusoidal, they are given 
by elliptic functions, but with a good accuracy one can write 0(t) ~ $sin(c\/A / $G?t), where 
c = ' ~ 0.85. The Universe at that time expands as at the radiation-dominated stage: 

ait) oc y/i. If one neglects the feedback of created 0-particles on the homogeneous field 4>(t), then 
its amplitude $(t) oc a _1 (t), so that a$ = const. Using a conformal time r], exact equation for 
quantum fluctuations S(p of the field can be reduced to the Lame equation. The results remain 
essentially the same if we use an approximate equation 



/fit 9t 
— ) (24) 
a{t) a(t) 

which leads to the Mathieu equation with A = ^^^2 + ^2 ~ ^^23,2 + 2.08, and q = ~ 1-04. 
Looking at the instability chart, we see that the resonance occurs in the second band, for k 2 ~ 
3Aa 2 $ 2 . The maximal value of the coefficient in this band for q ~ 1 approximately equals to 
0.07. As long as the backreaction of created particles is small, expansion of the Universe does not 
shift fluctuations away from the resonance band, and the number of produced particles grows as 
exp(2c AifeV / Aa$r/) ~ exp(^). 
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After the time interval ~ M p " 1 A~" 1 / 2 | In A|, backreaction of created particles becomes signifi- 
cant. The growth of the fluctuations (0 2 ) gives rise to a contribution 3A(0 2 ) to the effective mass 
squared of the field 0, both in the equation for (f)(t) and in Eq. fl2~4] ) for inhomogeneous modes. 
The stage of explosive reheating ends when (0 2 ) becomes greater than $ 2 . After that, $ 2 <C (0 2 ) 



and the effective frequency of oscillations is determined by the term y3A(0 2 ). The corresponding 

process is described by Eq. (4) with A{k) = 1 + 2q + ^Jb^sy, q = i^y- m this regime q <^ 1, 
and particle creation occurs in the narrow resonance regime in the second band with A ~ 4. 
Decay of the field in this regime is extremely slow: the amplitude $ decreases only by a factor 
t 1//12 faster that it would decrease without any decay, due to the expansion of the Universe only, 
i.e., $ oc t~ 7 / 12 P5| |. Reheating stops altogether when the presence of non-zero mass though 



still small as compared to y'3A(0 2 ) appears enough for the expansion of the Universe to drive a 
mode away from the narrow resonance. It happens when the amplitude $ drops up to a value 



In addition to this process, the field may decay to x-P ar tid es - This is the leading process 
for g 2 ^> A. The equation for Xk quanta has the same form as eq. (24) with the obvious change 
A — > g 2 /3. Initially parametric resonance is broad. The values of the parameter along the line 
A = 2q do not change monotonically, but typically for q ^> 1 they are 3 to 4 times greater than the 
parameter /i^ for the decay of the field <fi into its own quanta. Therefore, this pre-heating process 
is very efficient. It ends at the moment t ~ M^A" 1 / 2 ln(A/g 10 ) when $ 2 ~ (% 2 ) ~ g~ l \f\M 2 . 
The typical energy of created x-particles is E ~ (g 2 X) l l A M p . The following evolution is essentially 
the same as that described above for the case of a massive scalar field decaying into x-particles. 



Finally, let us consider the case with symmetry breaking. In the beginning, when the ampli- 
tude of oscillations is much greater than a, the theory of decay of the inflaton field is the same as 
in the case considered above. The most important part of pre-heating occurs at this stage. When 
the amplitude of the oscillations becomes smaller than m^/y/X and the field begins oscillating 
near the minimum of the effective potential at <fi = a, particle production due to the narrow 
parametric resonance typically becomes very weak. The main reason for this is related to the 
backreaction of particles created at the preceding stage of pre-heating on the rate of expansion of 
the Universe and on the shape of the effective potential ||43|| . However, importance of spontaneous 
symmetry breaking for the theory of reheating should not be underestimated, since it gives rise 
to the interaction term g 2 a(px 2 which is linear in cf). Such terms are necessary for a complete 
decay of the inflaton field in accordance with the perturbation theory (p|). 



Let us briefly summarize our results: 



1. In many models where decay of the inflaton field can occur in the purely bosonic sector 
the first stages of reheating occur due to parametric resonance. This process (pre-heating) is 
extremely efficient even if the corresponding coupling constants are very small. However, there is 
no explosive reheating in the models where decay of the inflaton field is necessarily accompanied 
by fermion production. 

2. The stage of explosive reheating due to a broad resonance typically is very short. Later the 



20 



resonance becomes narrow, and finally the stage of pre-heating finishes altogether. Interactions 
of particles produced at this stage, their decay into other particles and subsequent thermaliza- 
tion typically require much more time that the stage of pre-heating, since these processes are 
suppressed by the small values of coupling constants. 



3. The last stages of reheating typically can be described by the elementary theory of re- 



heating Pq| . However, this theory should be applied not to the decay of the original large and 
homogeneous oscillating inflaton field, but to the decay of particles produced at the stage of pre- 
heating, as well as to the decay of small remnants of the classical inflaton field. This makes a lot 
of difference, since typically coupling constants of interaction of the inflaton field with matter are 
extremely small, whereas coupling constants involved in the decay of other bosons can be much 
greater. As a result, the reheating temperature can be much higher than the typical temperature 
T r < 10 9 GeV which could be obtained neglecting the stage of parametric resonance [Q. On the 
other hand, such processes as baryon creation after inflation occur best of all outside the state of 
thermal equilibrium. Therefore, the stage of pre-heating may play an extremely important role 
in our cosmological scenario. Another consequence of the resonance effects is an almost instanta- 
neous change of equation of state from the vacuum-like one to the equation of state of relativistic 
matter. 



4. Reheating never completes in the theories where a single 0-particle cannot decay into other 
particles. This implies that reheating completes only if the theory contains interaction terms like 
4>ax 2 of (jnpip. In most cases the theories where reheating never completes contradict observational 
data. On the other hand, this result suggests an interesting possibility that the classical scalar 
fields (maybe even the inflaton field itself) may be responsible for the dark matter of the Universe 
even if they strongly interact with other matter fields. 



6 Topological defects as seeds for eternal inflation 



One of the main purposes of inflationary cosmology was to solve the primordial monopole problem. 
The solution was very simple: inflation exponentially increases the distance between monopoles 
and makes their density negligibly small. It was assumed that inflation neither can change internal 
properties of monopoles nor can it lead to their copious production. For example, in the first 
version of the new inflation scenario based on the £77(5) Coleman- Weinberg theory [BJ the Hubble 
constant during inflation was of the order 10 10 GeV, which is five orders of magnitude smaller 
than the mass of the X-boson Mx ~ 10 15 GeV. The size of a monopole given by ~ M^ 1 is five 
orders of magnitude smaller than the curvature of the Universe given by the size of the horizon 
if -1 . It seemed obvious that such monopoles simply could not know that the Universe is curved. 

This conclusion finds an independent confirmation in the calculation of the probability of 



spontaneous creation of monopoles during inflation. According to [49 , this probability is sup- 
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pressed by a factor of exp(— 2nm/H), where m is the monopole mass. In the model discussed 
above this factor is given by ~ 1CT 10 , which is negligibly small. 

Despite all these considerations, in the present section (see also ||28|| ) we will show that in many 
theories, including the SU(5) Coleman- Weinberg theory, monopoles, as well as other topological 
defects, do expand exponentially and can be copiously produced during inflation. The main 
reason is very simple. In the center of a monopole the scalar field <fi always vanishes. This means 
that monopoles always stay on the top of the effective potential at <ft = 0. When inflation begins, 
it makes the field <fi almost homogeneous and very close to zero near the center of the monopole. 
This provides excellent conditions for inflation inside the monopole. If inflation is supported 
by the vacuum energy V(0), the conditions for inflation remain satisfied inside the monopoles 
even after inflation finishes outside them. An outside observer will consider such monopoles as 
magnetically charged Reissner-Nordstrom black holes. A typical distance between such objects 
will be exponentially large, so they will not cause any cosmological problems. However, each such 
magnetically charged black hole will have an inflationary universe inside it. And, as we are going 
to argue, each such inflationary universe will contain many other inflating monopoles, which after 
inflation will look like black holes with inflationary heart. We will call such configurations fractal 
monopoles, or, more generally, fractal topological defects. 

To explain the basic idea in a more detailed way, we will consider first the simple model with 
the Lagrangian 

L = \{d,<P?-\^-^)\ (25) 

where is a real scalar field. Symmetry breaking in this model leads to formation of domains with 
(f) = ±7], where rj — These domains are divided by domain walls (kinks) which interpolate 
between the two minima. Let us for a moment neglect gravitational effects. In this case one can 
easily obtain a solution for a static domain wall in the yz plane: 

4> = i] tanh 




Now let us see whether our neglect of gravitational effects is reasonable. The potential energy 
density in the center of the kink (|26|) at x = is equal to ^ 4 , the gradient energy is also 
equal to \rf. This energy density remains almost constant at \x\ <C mr x = and then it 
rapidly decreases. Gravitational effects can be neglected if the Schwarzschild radius r g = 2M 
corresponding to the distribution of matter with energy density p = \rf and radius R ~ mr 1 is 
much smaller than R. Here M = ^-pR 3 . This condition suggests that gravitational effects can 
be neglected for r\ <C j-. In the opposite case, 



V Z 7T- , (27) 

Z7T 



gravitational effects can be very important. 



This result remains true for other topological defects as well. For example, recently it was 
shown that monopoles in the theory with the scale of spontaneous symmetry breaking t] > form 
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Reissner-Nordstrom black holes pD| . It remained unnoticed, however, that the same condition, 
r] > 1, is simultaneously a condition of inflation at <C r\ in the model (|25|) . 



Indeed, inflation occurs at <C rj in the model fl25| ) if the curvature of the effective potential 
V{4>) at <C rj is much smaller than 3iJ 2 , where = \f^V 2 is the Hubble constant supported 
by the effective potential ||. This gives m 2 <C 27rAr/ 4 , which leads to the condition almost exactly 
coinciding with ( |T5| ) : i] ^> / \^2~k. 

This coincidence by itself does not mean that domain walls and monopoles in the theories 
with 77 ^> 1 / \[2~n will inflate. Indeed, inflation occurs only if the energy density is dominated 
by the vacuum energy. As we have seen, for the wall (|26|) this was not the case: gradient energy 
density for the solution (^) near x = is equal to the potential energy density. However, this is 
correct only after inflation and only if gravitational effects are not taken into account. 

At the initial stages of inflation the field is equal to zero. Even if originally there were 
any gradients of this field, they rapidly become exponentially small. Each time At = H -1 new 
perturbations with the amplitude H/\/2tt and the wavelength ~ H^ 1 are produced, but their 
gradient energy density ~ H 4 is always much smaller than V(4>) for V(<fi) <C 1 P^l . 

Formation of domain walls in this scenario can be explained as follows. Assume that in the 
beginning the field inside some domain of initial size O^H^ 1 ) was equal to zero. It did not 
move classically at that time, since V'(0) = 0. However, within the time At = H~ l fluctuations 
of the field were generated, which looked like sinusoidal waves with an amplitude H/ \/2n 
and wavelength ~ H~ x . During this time the original domain grows in size e times, its volume 
grows e 3 ~ 20 times. Therefore it becomes divided into 20 domains of a size of the horizon 
H~ x . Evolution of the field inside each of them occurs independently of the processes in the 
other domains (no-hair theorem for de Sitter space). In a half of these domains scalar field will 
have average value ~ +H/2n, in other domains it will have average value ~ —H/2ir. After 
that, the field begins its classical motion. In those domains where the resulting amplitude of 
the fluctuations is positive, the field moves towards the minimum of the effective potential at 
= +rj] in other domains it moves towards = — rj. These domains become separated by the 
domain walls with = 0. However, as we already mentioned, at this stage the gradient energy 
near these walls is similar to H A ~ V 2 (<p), which is much smaller than V(<p) for V(<p) 1. Thus, 
domain walls continue expanding exponentially in all directions. 

Domain wall formation is not completed at this point yet. Indeed, the field may jump back 
from the region with > to < 0. However, the probability of such jumps is smaller than 1/2, 
and it becomes even much smaller with a further growth of |0|. What is more important, the 
typical wavelength of new fluctuations remains H" 1 , whereas the previously formed domains with 
~ ±H/2n continue growing exponentially. Therefore the domains with negative produced by 
the jumps of the field will appear inside domains with positive 0, but they will look like small 
islands with < inside the sea with > 0. This means that the subsequent stages of inflation 
cannot destroy the originally produced domain walls; they can only produce new domain walls 
on a smaller length scale. These new walls will be formed only in those places where the scalar 
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field is sufficiently small for the jumps with the change of the sign of the field to be possible. 
Therefore the new walls will be created predominantly near the old ones (where = 0), thus 
forming a fractal domain wall structure. 



Similar effects occur in more complicated models where instead of a discrete symmetry — > — 
we have a continuous symmetry. For example, instead of the model (p5|) one can consider a model 

L = d^*d^-xU*<f>-^\ , (28) 

where is a complex scalar field, = 4= (0i + i(f>2)- Spontaneous breaking of the U(l) symmetry 
in this theory may produce global cosmic strings. Each string contains a line with = 0. Outside 
this line the absolute value of the field increases and asymptotically approaches the limiting 
value ^0^ + 02 = 77. This string will be topologically stable if the isotopic vector (0!(x), 2 (x)) 
rotates by 2rm when the point x takes a closed path around the string. 



The next step is to consider a theory with 0(3) symmetry, 

L = \{d^f-\^-i?)\ (29) 

where is a vector (01,02,03). This theory admits global monopole solutions. The simplest 
monopole configuration contains a point x = with 0(0) = surrounded by the scalar field 
0(x) oc x. Asymptotically this field approaches regime with 2 (x) = rf. 



The most important feature of strings and monopoles is the existence of the points where 
= 0. Effective potential has an extremum at = 0, and if the curvature of the effective potential 
is smaller than H 2 = ®ZLKM ; space around the points with = will expand exponentially, just 
as in the domain wall case considered above. 



Now we will add gauge fields. We begin with the Higgs model, which is a direct generalization 
of the model @: 

1 

T 




L = D^cjf D^4> - -F^F^ - A ( 0*0 - ] . (30) 



Here is a covariant derivative of the scalar field, which in this simple case (Abelian theory) 
is given by — ieA^. In this model strings of the scalar field contain magnetic flux $ = 27r/e. 
This flux is localized near the center of the string with <j)(x) = 0, for the reason that the vector 



field becomes heavy at large 0, see e.g. ||14|| . However, if inflation takes place inside the string, 
then the field becomes vanishingly small not only at the central line with 0(x) = 0, but even 
exponentially far away from it. In such a situation the flux of magnetic field will not be confined 
near the center of the string. The thickness of the flux will grow together with the growth of the 
Universe. Since the total flux of magnetic field inside the string is conserved, its strength will 
decrease exponentially, and very soon its effect on the string expansion will become negligibly 
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small. A deep underlying reason for this behavior is the conformal invariance of massless vector 
fields. Energy density of such fields decreases as a(t)~ 4 , where a(t) is a scale factor of the Universe. 
The final conclusion is that the vector fields do not prevent inflation of strings. The condition 
for inflation to occur inside the string remains the same as before: the curvature of the effective 
potential at = should be smaller than H 2 = 87r ^°- > . 

The final step is to consider magnetic monopoles. A simplest example is given by the 0(3) 
theory 

L = \\(Dj\ 2 - \f; u F^ - - rff . (31) 

Global monopoles of the theory ( p8[ ) become magnetic monopoles in the theory fl3"0|). They also 
have = in the center. Vector fields in the center of the monopole are massless (g<p = 0). 
During inflation these fields exponentially decrease, and therefore they do not affect inflation of 
the monopoles. 

We should emphasize that even though the field (f) around the monopole during inflation is 
very small, its topological charge is well defined, it cannot change and it cannot annihilate with 
the charge of other monopoles as soon as the radius of the monopole becomes greater than H . 
However, an opposite process is possible. Just as domain walls can be easily produced by quantum 
fluctuations near other inflating domain walls, pairs of monopoles can be produced in the vicinity 
of an inflationary monopole. The distance between these monopoles grow exponentially, but the 
new monopoles will appear in the vicinity of each of them. This picture has been confirmed by 
the results of computer simulation of this process . 

Note that in the simple models discussed above inflation of monopoles occurs only if sponta- 
neous symmetry breaking is extremely strong, r\ > 1. However, this is not a necessary condition. 
Our arguments remain valid for all models where the curvature of the effective potential near 
4> = is smaller than the Hubble constant supported by ^(0). This condition is satisfied by 
all models which were originally proposed for the realization of the new inflationary universe 
scenario. In particular, the monopoles in the SU(5) Coleman- Weinberg theory also should ex- 
pand exponentially. However, to make sure that we do not miss something important, we should 
resolve the paradox formulated in the very beginning of this section: The Hubble constant H 
during inflation in the SU(5) Coleman- Weinberg theory is much smaller than the mass of the 
vector field Mx, which is usually related to the size of the monopole. In such a situation inflation 
of the interior of the monopole does not seem possible. 

This objection is very similar to the argument that the domain wall should not inflate since 
its gradient energy density is equal to its potential energy density. The resolution of the paradox 
is very similar. The effective mass of the vector field Mx ~ grj ~ 10 15 GeV can determine an 
effective size of the monopole only after inflation. Effective mass of the vector field Mx{4>) ~ g<p 
is always equal to zero in the center of the monopole. Once inflation begins in a domain of a 
size OlH^ 1 ) around the center of the monopole, it expels vectors fields away from the center and 
does not allow them to penetrate back as far as inflation continues. 
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As we already discussed, it in many versions of inflationary theory some parts of the Universe 
expand without end. Our results reveal a very nontrivial role which topological defects may play 
in this process. For example, in the theory with breaking of a discrete symmetry — > — the 
Universe in the source of its evolution becomes divided into many thermalized regions divided by 
exponentially expanding domain walls unceasingly producing new inflating walls. It is important 
that inflation never stops near the walls. Some of these walls could collapse, eating the island of 
thermalized phase inside them. However, it does not seem possible. Due to the no-hair theorem 
for de Sitter space, the part of the Universe near the wall (i.e. near = 0) lives by its own laws 
and continue expanding exponentially all the time. Thus, the walls become indestructible sources 
of eternal inflation, which produce new inflating walls, which produce new walls, etc. Far away 
from the walls, the field relaxes near = ±77. However, due to the exponential expansion of 
space near the walls, they always remain exponentially thick and never approach the thin wall 
solution described by eq. (|13|). 



The picture of eternally inflating Universe consisting of islands of thermalized phase sur- 
rounded by inflating domain walls is very similar to a picture which appears in old inflation and 
in the versions of new inflation where the state = is sufficiently stable. However, in our model 
the main reason for this behavior is purely topological ||28|| . If inflation in this theory begins at 
large 0, as in the simplest models of chaotic inflation, then the Universe will not consist of islands 
of thermalized phase surrounded by de Sitter space. On the contrary, it will consist of islands 
of inflationary universe surrounded by the thermalized phase ||32|| . Nevertheless, inflation in this 



case will go eternally due to the process of self-reproduction of inflationary domains with large 
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The results of computer simulations performed in |28| suggest that in the theories with con- 
tinuous symmetry breaking, which allow existence of inflating strings and monopoles, the global 
structure of the Universe is similar to the structure of the Universe in the simplest versions of 
chaotic inflation. The Universe will consist of islands of inflationary phase associated with in- 
flating monopoles (or of lines of inflationary phase associated with inflating strings) surrounded 
by thermalized phase. Since monopoles are topologically stable, and field is equal to zero in 
their centers, they form indestructible sources of inflation. The structure of space-time near each 



such monopole is very complicated; it should be studied by the methods developed in ]5I] for 



description of a bubble of de Sitter space immersed into vacuum with vanishing energy density. 
Depending on initial conditions, many possible configurations may appear. The simplest one is 
the monopole which looks like a small magnetically charged black hole from the outside, but 
which contains a part of exponentially expanding space inside it. This is a wormhole configura- 
tion similar to those studied in p2"| — [p4] . The physical interpretation of this configuration can 



be given as follows [51, 54]. An external observer will see a small magnetically charged Reissner- 



Nordstrom black hole. This part of the picture is consistent with the results obtained in [5T5]. On 



the other hand, an observer near the monopole will see himself inside an eternally inflating part 
of the Universe. These two parts of the Universe will be connected by a wormhole. Note that the 
wormhole connecting inflationary magnetic monopole to our space cannot evaporate unless the 
inflationary universe can loose its magnetic charge. 



At the quantum level the situation becomes even more interesting. As we already mentioned, 
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fluctuations of the field near the center of a monopole are strong enough to create new regions 
of space with = 0, some of which will become monopoles. After a while, the distance between 
these monopoles becomes exponentially large, so that they cannot annihilate. This process of 
monopole-antimonopole pair creation produces a fractal structure consisting of monopoles created 
in the vicinity of other monopoles. Each of these monopoles in the process of its further evolution 
will evolve into a black hole containing inflationary universe containing many other monopoles, 
etc. That is why we call them fractal monopoles. 



Note that one of the main motivations for the development of inflationary cosmology was a 
desire to get rid of primordial monopoles. However, the original version of the new inflation- 
ary universe scenario based on the theory of high temperature phase transitions did not work 
particularly well, since it required rather unnatural initial conditions ||. Fortunately, the same 
potentials which gave rise to new inflation can be used in the context of the chaotic inflation 
scenario, which does not require the whole Universe being in a state of thermal equilibrium with 
= 0. If these initial conditions can be realized at least in one domain of initial size 0(i/ _1 ), the 
Universe enters the stage of eternal inflation. Of course, this possibility is still rather problematic 
since in the models used in the new inflationary universe scenario the size 0(if _1 ) is many orders 
of magnitude greater than the Planck length. For example, the probability of quantum creation of 
inflationary universe in the new inflationary universe scenario is suppressed by an infinitesimally 
small factor ex p(~ 8 y 3 (o) ) ~ 10~ 10l ° W^\ - Therefore it is much easier for inflation to begin and con- 
tinue eternally in the simplest chaotic inflation models where inflation is possible even at H ~ 1, 
V{4>) ~ 1 |S2|, [IIJ. Still the problem of initial conditions in the models with the potentials used 
in the new inflationary scenario does not look that bad if one takes into account that inflation in 
these models also goes on without end [3D], ^TJ, and eventually we become generously rewarded for 
our initial problems by the infinite growth of the total volume of inflationary domains. Moreover, 
as it was shown in p2|| , the problem of initial conditions for the models with the potentials used in 
the new inflationary scenario can be easily solved by a preceding stage of chaotic inflation. Now 
we encounter a new amazing twist of this scenario. Eternal inflation begins if initially we have at 
least one inflating monopole. Such monopoles exponentially expand, create new space suitable 
for life of our type to appear there, and then these monopoles effectively expel themselves from 
the space where inflation ended by continuous creation of exponentially expanding space around 
them. The possibility that monopoles by themselves can solve the monopole problem is so simple 
that it certainly deserves further investigation. 



7 Stationary Universe 



Even though eternal inflation due to topological defects looks very interesting, it is much less 
general than the eternal inflation which occur because of quantum fluctuations in many other 



theories including the simplest models of chaotic inflation with the potentials ! y 1 2 or j0 4 ED 



In fact, the effect of self- reproduction of inflationary universe occurs almost in all versions of in- 



27 



flationary cosmology (even though some exceptions are possible ||34|| ). This forces us to reconsider 
many features of the big bang cosmology. 

The first models of inflation were based on the standard assumption of the big bang theory 
that the Universe was created at a single moment of time in a state with the Planck density, and 
that it was hot and large (much larger than the Planck scale M~ x = 1) from the very beginning. 
The success of inflation in solving internal problems of the big bang theory apparently removed 
the last doubts concerning the big bang cosmology. It remained almost unnoticed that during the 
last ten years the inflationary theory has broken the umbilical cord connecting it with the old big 
bang theory, and acquired an independent life of its own. For the practical purposes of description 
of the observable part of our Universe one may still speak about the big bang. However, if one 
tries to understand the beginning of the Universe, or its end, or its global structure, then some 
of the notions of the big bang theory become inadequate. 

For example, already in the first version of the chaotic inflation scenario there was no 
need to assume that the whole Universe appeared from nothing at a single moment of time 
associated with the big bang, that the Universe was hot from the very beginning and that the 
inflaton scalar field <fi which drives inflation originally occupied the minimum of its potential 
energy. Later it was found that if the Universe contains at least one inflationary domain of a 
size of horizon ( L h- region') with a sufficiently large and homogeneous scalar field <p, then this 
domain will permanently produce new /i-regions of a similar type. In other words, instead of a 
single big bang producing a one-bubble Universe, we are speaking now about inflationary bubbles 
producing new bubbles, producing new bubbles, ad infinitum. In this sense, inflation is not a short 
intermediate stage of duration ~ 10~ 35 seconds, but a self-regenerating process, which occurs in 
some parts of the Universe even now, and which will continue without end. 

Thus, recent developments of inflationary theory have considerably modified our cosmological 
paradigm ||. Now we must learn how to formulate physical questions in the new context. For 
example, in a homogeneous part of the Universe there is a simple relation between the density of 
matter and time. However, on a very large scale the Universe becomes extremely inhomogeneous. 
Its density, at the same 'cosmic time', varies anywhere from zero to the Planck density. Therefore 
the question about the density of the Universe at the time 10 10 years may not have any definite 
answer. Instead of addressing such questions we should study the distribution of probability of 
finding a part of the Universe with given properties, and find possible correlations between these 
properties. 

It is extremely complicated to describe an inhomogeneous Universe and to find the corre- 
sponding probability distribution. Fortunately, there exists a particular kind of stationarity of 
the process of the Universe self-reproduction which makes things more regular. Due to the no- 
hair theorem for de Sitter space, the process of production of new inflationary domains occurs 
independently of any processes outside the horizon. This process depends only on the values of 
the fields inside each /i-region of radius H~ l . Each time a new inflationary /i-region is created 
during the Universe expansion, the physical processes inside this region will depend only on the 
properties of the fields inside it, but not on the 'cosmic time' at which it was created. 
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In addition to this most profound stationarity, which we will call local stationarity there may 
also exist some simple stationary probability distributions which may allow us to say, for example, 
what the probability is of finding a given field at a given point. To examine this possibility one 
should consider the probability distribution P c (<p,t\x), which describes the probability of finding 
the field at a given point at a time t, under the condition that at the time t = the field at 
this point was equal to x- The same function may also describe the probability that the scalar 
field which at time t was equal to 0, at some earlier time t = was equal to x- 



The probability distribution P c has been studied by many authors, see e.g. refs. |T0|, |55], j56f . 
Our investigation of this question has shown that in all realistic inflationary models the proba- 
bility distribution P c (0, t\x) is not stationary [ffD . The reason is very simple. The probability 



distribution P c is in fact the probability distribution per unit volume in comoving coordinates 
(hence the index c in P c ), which do not change during expansion of the Universe. By considering 
this probability distribution we neglect the main source of the self-reproduction of inflationary 
domains, which is the exponential growth of their volume. Therefore, in addition to P c , we in- 
troduced the probability distribution P p (<f), t\x), which describes the probability to find a given 
field configuration in a unit physical volume [0. In this section we will show that under certain 



conditions the stationary probability distribution P p ((f),t\x) does exist P2[ . In such situations we 
will speak about global stationarity. 



First of all we should remember some details of stochastic approach to inflation. Let us con- 
sider the simplest model of chaotic inflation based on the theory of a scalar field minimally 
coupled to gravity, with the effective potential V(4>). If the classical field is sufficiently homo- 
geneous in some domain of the Universe, then its behavior inside this domain is governed by the 
equation 3H<p = -dV/d<p } where H 2 = W*. 

Investigation of these equations has shown that in all power-law potentials V(<f)) ~ n inflation 
occurs at > e ~ n/6. In the theory with an exponential potential V(<ft) ~ e°"^ inflation ends 
only if we bend down the potential at some point e ; for definiteness we will take e = in this 
theory. 

Inflation stretches all initial inhomogeneities. Therefore, if the evolution of the Universe were 
governed solely by classical equations of motion, we would end up with an extremely smooth 
Universe with no primordial fluctuations to initiate the growth of galaxies. Fortunately, new 
density perturbations are generated during inflation due to quantum effects. The wavelengths of 
all vacuum fluctuations of the scalar field grow exponentially in the expanding Universe. When 
the wavelength of any particular fluctuation becomes greater than H~ x , this fluctuation stops 
oscillating, and its amplitude freezes at some nonzero value <50(x) because of the large friction 
term 3if in the equation of motion of the field 0. The amplitude of this fluctuation then remains 
almost unchanged for a very long time, whereas its wavelength grows exponentially. Therefore, 
the appearance of such a frozen fluctuation is equivalent to the appearance of a classical field 
S(f)(x) that does not vanish after averaging over macroscopic intervals of space and time. 

Because the vacuum contains fluctuations of all wavelengths, inflation leads to the creation 
of more and more perturbations of the classical field with wavelengths greater than if -1 . The 
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average amplitude of such perturbations generated during a time interval H 1 (in which the 
Universe expands by a factor of e) is given by 

1^)1^- (32) 

The phases of each wave are random. Therefore, the sum of all waves at a given point fluctuates 
and experiences Brownian jumps in all directions in the space of fields. 

The standard way to describe the stochastic behavior of the inflaton field during the slow- 
rolling stage is to coarse-grain it over /i-regions and consider the effective equation of motion of 
the long- wavelength field [ 55fl : 



Here H = ^8nV/3 , £ (t) is the effective white noise generated by quantum fluctuations, which 
leads to the Brownian motion of the classical field 4>. 

This Langevin equation leads to two stochastic equations for the probability distribution 
P c ((j),t\x)- The first one is called the backward Kolmogorov equation, 

—dT~ = a* [-*r tlx) )~ mdK m ' tlx) • (34) 

In this equation one considers the value of the field <ft at the time i as a constant, and finds the 
time dependence of the probability that this value was reached during the time t as a result of 
diffusion of the scalar field from different possible initial values \ = (f)(0). 

The second equation is the adjoint to the first one; it is called the forward Kolmogorov 
equation, or the Fokker-Planck equation flSTj . 



One may try to find a stationary solution of equations fl34|), (p5|), assuming that 9Pc ^ x ^ = 0. 
The simplest stationary solution (subexponential factors being omitted) would be 

PM - t]x) ~^(^)- e * p (-mxi) ■ (36) 

This function is extremely interesting. Indeed, the first term in fl36|) is equal to the square of the 



Hartle-Hawking wave function of the Universe |57|], whereas the second one gives the square of 



the tunneling wave function [ 373 1 ! 



At first glance, this result gives a direct confirmation and a simple physical interpretation 
of both the Hartle-Hawking wave function of the Universe and the tunneling wave function. 
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However, in all realistic cosmological theories, in which V{4>) = at its minimum, the Hartle- 
Hawking distribution exp ( g^j j is not normalizable. The source of this difficulty can be easily 
understood: any stationary distribution may exist only due to compensation of the classical flow 
of the field downwards to the minimum of V(0) by the diffusion motion upwards. However, 
diffusion of the field discussed above exists only during inflation. Thus, there is no diffusion 
motion upwards from the region < <p e . Therefore all solutions of equation (|35|) with the proper 
boundary conditions at = e (i.e. at the end of inflation) are non-stationary (decaying) [ |iO|j . 

The situation with the probability distribution P p is much more interesting and complicated. 
As was shown in |]T(J its behavior depends strongly on initial conditions. If the distribution P p 



was initially concentrated at < 0*, where 0* is some critical value of the field, then it moves 
towards small in the same way as P c , i.e. it cannot become stationary. On the other hand, if 
the initial value of the field is larger than 0*, the distribution moves towards larger and larger 
values of the field 0*, until it reaches the field P , at which the effective potential of the field 
becomes of the order of Planck density M p = 1, where the standard methods of quantum field 
theory in a curved classical space are no longer valid. 



Some further steps towards the solution of this problem were made by Nambu and Sasaki [58 



and Mijic 59 . Their papers contain many important results and insights. However, Mijic E| 



did not have a purpose to obtain a complete expression for the stationary distribution P p (0, t\x)- 
The corresponding expressions were obtained for various types of potentials V(4>) in [58R. Unfor- 



tunately, according to the stationary distribution P p ((j),t\x) is almost entirely concentrated 



at ^> p , i.e. at V((f)) ^> 1, where the methods used in |58| are inapplicable. 



We will continue this investigation by writing the system of stochastic equations for P p . These 
equations can be obtained from eqs. (0), ( p5|) by adding the term 3HP P , which appears due 
to the growth of physical volume of the Universe by the factor 1 + 3iJ(0) dt during each time 



interval dt 0, [H], [3 

dP p _ 1 H 3 / 2 ( X ) d (H 3 / 2 ( X )dP p \ V'{ X )dP 1 



v 



dt 2 2tt d X V 2vr d X J SH( X ) d X 



+ 3H( X )P P , (37) 



dt 2d(p\ 2tt d<p V 2tt " v > 3H(<j>) v , 

To find solutions of these equations one must specify boundary conditions. Behavior of solutions 
typically is not very sensitive to the boundary conditions at e ; it is sufficient to assume that the 



diffusion coefficient (and, correspondingly, the first terms in the r.h.s. of equations (37), (|38|) ) 
vanish for < e . 

The conditions at the Planck boundary = P play a more important role. Investigation of 
these conditions poses many difficult problems. First of all, our diffusion equations are based on 
the semiclassical approach to quantum gravity, which breaks down at super-Planckian densities. 
Secondly, the shape of the effective potential may be strongly modified by quantum effects at 
~ (p p . Finally, our standard interpretation of the probability distribution P c and P p breaks 
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down at V{4>) > 1, since at super-Planckian densities the notion of a classical scalar field in 
classical space-time does not make much sense. 

However, these problems by themselves suggest a possible answer. Inflation happens only in 
theories with very flat effective potentials. At the Planck density nothing can protect the effective 
potential from becoming steep. Hence, one may expect that inflation ceases to exist at <ft > (ftp, 
which leads to the boundary condition 

P p (<f> p ,t\x) = P p (<j>,t\ Xp )=0, (39) 

where V(<P P ) = V( Xp ) = 0(1). 

There exist some other reasons discussed in [^, why inflation and/or self- reproduction 

of inflationary domains may cease to exist at super-Planckian densities. We will just assume that 
this is the case and impose the boundary conditions fl39|) on P p . 

One may try to obtain solutions of equations (0), ( p8|) in the form of the following series of 
biorthonormal system of eigenfunctions of the pair of adjoint linear operators (defined by the left 
hand sides of the equations below): 

oo 

^(<Mx) = E eAs >s(xK(</>) • (40) 



s=l 



Indeed, this gives us a solution of eq. (^) if 

1 # 3/2 (x) d ( H*I*{ X ) d , , A V'{x) d 

- 2 —^q- x \-^Tq^M ) - 3H[xjdx + ' MX) = K MX) ■ (41) 

and 



The orthonormality condition reads 



^s{x)^j(x)dX = $sj • (43) 



In our case (with regular boundary conditions) one can easily show that the spectrum of Xj is 
discrete and bounded from above. Therefore the asymptotic solution for P p (<f>,t\x) (in the limit 
t — > oo) is given by 

PpOMM =e Xlt Mx)*i(<!>) ■ (l + 0(e-^- x ^)) . (44) 
Here V'i(x) ls the only positive eigenfunction of eq. (flTl), Ai is the corresponding (real) eigenvalue, 



and 7Ti(0) is the eigenfunction of the conjugate operator ( fi2]) with the same eigenvalue \\. Note, 
that Ai is the highest eigenvalue, Re (Ai — A 2 ) > 0. This is the reason why the asymptotic 
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equation (H) is valid at large t. We have found PS] that in realistic theories of inflation a typical 
time of relaxing to the asymptotic regime, At ~ (Ai — A2) -1 , is extremely small. Typically it is 
only about a few thousands Planck times, i.e. about 10 -40 sec. This means that the normalized 
distribution 

P p (<M| X )=e- Al *P p (Mx) ( 45 ) 
rapidly converges to the time-independent normalized distribution 

PMx) = Pp(<j>,t^ oo\ X ) = MX) ?ri(0) ■ (46) 

It is this stationary distribution that we were looking for. Because the growing factor e~ Xlt is the 
same for all (and x)> one can use Pp instead of P p for calculation of all relative probabilities. 
In particular, P p (<f>\x) gives us the fraction of the volume of the Universe occupied by the field 
0, under the condition that the corresponding part of the Universe at some time in the past 
contained the field \. The remaining problem is to find the functions ipi(x) and 7Ti(0), and to 
check that all assumptions about the boundary conditions which we made on the way to eq. ([44"D 
are actually satisfied. 

We have solved this problem for chaotic inflation in a wide class of theories including the 
theories with polynomial and exponential effective potentials V{4>) and found the corresponding 
stationary distributions f52 |. Here we will present some of our results for the theories ^0 4 and 

V n e a<i >. 



Solution of equations ( J4T| ) and ( f4"2"D for ipx(x) and 7Ti(0) in the theory |0 4 shows that these 
functions are extremely small at ~ 0e and X ~ Xe- They grow at large and Xi then rapidly 
decrease, and vanish at = x — 4>p- With a decrease of A the solutions become more and 
more sharply peaked near the Planck boundary. (The functions ipi and 7Ti for the exponential 
potential have a similar behavior, but they are less sharply peaked near <f> p .) A detailed discussion 
of these solutions will be contained in [p2| . The eigenvalues Ai corresponding to different coupling 



constants A are given by the following table: 



A 


1 


10- 1 


10- 2 


lO" 3 


10" 4 


10" 5 


10~ 6 


Ai 


2.813 


4.418 


5.543 


6.405 


7.057 


7.538 


7.885 



One can find also the second eigenvalue A2. For example, for A = 10~ 4 one gets A2 = 6.789. This 
means that for A = 10~ 4 the time of relaxation to the stationary distribution is At ~ (Ai — A2) _1 ~ 
4M" 1 ~ 10 -42 seconds — a very short time indeed. 

Note, that the distribution P p (<f),t\x) depends on and x ver y sharply. For example, one can 
show that in the theory ^0 4 

p p (<Mx)~Vi(x)<^ Ai • ( 47 ) 

for a sufficiently small (for < A -1 / 8 ) This is an extremely sharp dependence indeed. For 
example, for the realistic value A ~ 10 -13 one obtains P p {4>, r\x) ~ '0i(x)0 10 ! This result will 
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be important for us later, when we will be discussing predictions of quantum cosmology for the 
value of the cosmological parameter Q |35 . 



Note that the parameter Ai shows the speed of exponential expansion of the volume filled by a 
given field 4>. This speed does not depend on the field <fi, and has the same order of magnitude as the 



speed of expansion at the Planck density Indeed, Ai should be compared to 3H(<f>) = 2^6irV((j)) } 

which is equal to 2v / 6~7r »s 8.68 at the Planck density. It can be shown [j32j] that in the limit A — > 
the eigenvalue Ai also becomes equal to 2\/6tt ~ 8.681. The meaning of this result is very simple: 
in the limit A — > our solution becomes completely concentrated near the Planck boundary, and 
Ai becomes equal to 3if (</> p ). 

One should emphasize that the factor e Al * gives the rate of growth of the combined volume of all 
domains with a given field <p (o r of all domains containing matter with a given density) not only at 
very large <fi, where quantum fluctuations are large, but at small <ft as well, and even after inflation 
|3"2"|| . This result may seem absolutely unexpected, since the volume of each particular inflationary 
domain grows like e 3J? ^*, and after inflation the law of expansion becomes completely different. 
One should distinguish, however, between the growth of each particular domain, accompanied 
by a decrease of density inside it, and the growth of the total volume of all domains containing 
matter with a given (constant) density. In the standard big bang theory the second possibility 
did not exist, since the energy density was assumed to be the same in all parts of the Universe 
("cosmological principle"), and it was not constant in time. 

The reason why there is a universal expansion rate e Al * can be understood as follows. Because 
of the self-reproduction of the Universe there always exist many domains with <ft ~ <p p , and their 
combined volume grows almost as fast as e 3H ^ p ^. Then the field <p inside some of these domains 
decreases. The total volume of domains containing some small field grows not only due to 
expansion ~ g 3-ff W*, but mainly due to the unceasing process of expansion of domains with large 
4> and their subsequent rolling (or diffusion) towards small 0. 

But what about the field which is already at the Planck boundary? Why do the correspond- 
ing domains not grow exactly with the Planckian Hubble constant H(<j) p ) = 2\/67r/3? It happens 
partially due to diffusion and slow rolling of the field towards smaller <p. However, the leading 
effect is the destructive diffusion towards the space-time foam with <fi > <p p . One may visualize 
this process by painting white all domains with V(4>) < 1, and by painting black domains filled 
by space-time foam with V{4>) > 1. Then each time if _1 (0 p ) the volume of white domains with 
(f) ~ cj) p grows approximately e 3 times, but some 'black holes' appear in these domains, and, as a 
result, the total volume of white domains increases only e 3Al / 2v/ ^ times. This suggests (by anal- 
ogy with [62]) calling the factor dj = 3Ai/2\/67r 'the fractal dimension of classical space-time', 



or 'the fractal dimension of the inflationary Universe'. However, one should keep in mind that 
the fractal structure of the inflationary Universe in the chaotic inflation scenario in general is 
more complicated than in the new or old inflation and cannot be completely specified just by one 



fractal dimension 32 



The distribution P p ((p\x) = i>i{x) ^litf") which we have obtained does not depend on time t. 
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However, in general relativity one may use many different time parametrizations, and the same 
physics can be described differently in different 'times'. One of the most natural choices of time 
in the context of stochastic approach to inflation is the time r = In °^o) = / H((f)(x,t),t) dt 
55], [56j . Here a (x, t) is a local value of the scale factor in the inflationary Universe. By using this 



time variable, we were able to obtain not only numerical solutions to the stochastic equations, 
but also simple asymptotic expressions describing these solutions. For example, for the theory 
j(p 4 both the eigenvalue Ai and the 'fractal dimension' df (which in this case refers both to the 
Planck boundary at (j) p and to the end of inflation at e ) are given by = Ai ~ 3 — 1.1 \/X, and 
the stationary distribution is 

3 11 

p„(<Mx) ~ ex p(-M7^y) { v(x) + o.4 ~ TV • ^ exp (~' r (3 ~ Al¥ ') 

~ exp (-2^i) (a^tt^ - n) ■ * M-^ 2 ) ■ («) 

Note that the first factor coincides with the square of the tunneling wave function []nj! This 
expression is valid in the whole interval from e to <p p and it correctly describes asymptotic 
behavior of P p (4>, t\x) both at x ~ Xe and at x ~ Xv 

A similar investigation can be carried out for the theory V(<fi) = V Q e a9 \ The corresponding 
solution is 

W,T|X) ~ exp(-^) ( ^ - 1) • - 1) V*«) . (49) 

This expression gives a rather good approximation for P p (4>,t\x) for all <j) and x- 



The main result of our work is that under certain conditions the properties of our Universe 
can be described by a time-independent probability distribution, which we have found for theories 
with polynomial and exponential effective potentials. A lot of work still has to be done to verify 
this conclusion, see |52J. However, once this result is taken seriously, one should consider its 
interpretation and rather unusual implications. 

When making cosmological observations, we study our part of the Universe and find that in 
this part inflation ended about t e ~ 10 10 years ago. The standard assumption of the first models 
of inflation was that the total duration of the inflationary stage was At ~ 10~ 35 seconds. Thus 
one could come to an obvious conclusion that our part of the Universe was created in the big 
bang, at the time t e + At ~ 10 10 years ago. However, in our scenario the answer is quite different. 

Let us consider an inflationary domain which gave rise to the process of self-reproduction 
of new inflationary domains. For illustrative purposes, one can visualize self-reproduction of 
inflationary domains as a branching process, which gives a qualitatively correct description of 
the actual physical process we consider. During this process, the first inflationary domain of 
initial radius ~ within the time splits into e 3 ~ 20 independent inflationary 

domains of similar size. Each of them contains a slightly different field 0, modified both by 
classical motion down to the minimum of V(0) and by long-wavelength quantum fluctuations of 
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amplitude ~ H/2it. After the next time step if _1 (0), which will be slightly different for each of 
these domains, they split again, etc. The whole process now looks like a branching tree growing 
from the first (root) domain. The radius of each branch is given by H~ x ] the total volume of 
all domains at any given time t corresponds to the 'cross-section' of all branches of the tree at 
that time, and is proportional to the number of branches. This volume rapidly grows, but when 
calculating it, one should take into account that those branches, in which the field becomes larger 
than (j) p , die and fall down from the tree, and each branch in which the field becomes smaller than 
4> e , ends on an apple (a part of the Universe where inflation ended and life became possible). 

One of our results is that even after we discard at each given moment the dead branches 
and the branches ended with apples, the total volume of live (inflationary) domains will continue 
growing exponentially, as e Al *. What is even more interesting, we have found that very soon 
the portion of branches with given properties (with given values of scalar fields, etc.) becomes 
time-independent. Thus, by observing any finite part of a tree at any given time t one cannot 
tell how old the tree is. 



To give a most dramatic representation of our conclusions, let us see where most of the apples 
grow. This can be done simply by integrating e Xlt from t = to t = T and taking the limit as 
T — > oo. The result obviously diverges at large T as Aj" 1 e AlT , which means that most apples 
grow at an indefinitely large distance from the root. In other words, if we ask what is the total 
duration of inflation which produced a typical apple, the answer is that it is indefinitely large. 



This conclusion may seem very strange. Indeed, if one takes a typical point in the root domain, 
one can show that inflation at this point ends within a finite time At ~ 10 -35 seconds. This is 
a correct (though model-dependent) result which can be confirmed by stochastic methods, using 
the distribution P c ((p,At\x) HlOU - How could it happen that the duration of inflation was any 
longer than 10~ 35 seconds? 

The answer is related to the choice between P c and P p , or between roots and fruits. Typical 
points in the root domain drop out from the process of inflation within 10~ 35 seconds. The number 
of those points which drop out from inflation at a much later stage is exponentially suppressed, 
but they produce the main part of the total volume of the Universe. Note that the length of each 



particular branch continued back in time may well be finite |k| . However, there is no upper limit 
to the length of each branch, and, as we have seen, the longest branches produce almost all parts 
of the Universe with properties similar to the properties of the part where we live now. Since 
by local observations we can tell nothing about our distance in time from the root domain, our 
probabilistic arguments suggest that the root domain is, perhaps, indefinitely far away from us. 
Moreover, nothing in our part of the Universe depends on the distance from the root domain, 
and, consequently, on the distance from the big bang. 

Thus, inflation solves many problems of the big bang theory and ensures that this theory 
provides an excellent description of the local structure of the Universe. However, after making 
all kinds of improvements of this theory, we are now winding up with a model of a stationary 
Universe, in which the notion of the big bang loses its dominant position, being removed to the 
indefinite past. 
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8 Cosmological constant, Q < 1, and all that 



When inflationary theory was first formulated, we did not know how much it was going to influence 
our understanding of the structure of the Universe. We were happy that inflation provided an 
easy explanation of the homogeneity of the Universe. However, we did not know that the same 
theory simultaneously predicts that on the extremely large scale the Universe becomes entirely 
inhomogeneous, and that this inhomogeneity is good, since it is one of the manifestations of the 
process of self-reproduction of inflationary universe. 

The new picture of the Universe which emerges now is very unusual, and we are still in the 
process of learning how to ask proper questions in the context of the new cosmological paradigm. 
Previously we assumed that we live in a Universe which has the same properties everywhere 
( "cosmological principle" ). Then one could make a guess about the most natural initial conditions 
in the Universe, and all the rest followed almost automatically. Now we learned that even if one 
begins with a non- uniform Universe, later it becomes extremely homogeneous on a very large 
scale. However, simultaneously it becomes absolutely non-uniform on a much greater scale. The 
Universe becomes divided into different exponentially large regions where the laws of low-energy 
physics become different. In certain cases the relative fraction of volume of the Universe in a 
state with given fields or with a given density does not depend on time, whereas the total volume 
of all parts of the Universe continues growing exponentially. 

In this situation the concept of naturalness becomes much more ambiguous than before, the 
problem of introduction of a proper measure of probability becomes most important. One of 
the most natural choices of such measure is given by the probability distribution P p (4>, t\x)- The 
hypothesis behind this proposal is that we are typical, and therefore we live in those parts of the 
Universe where most other people do. The total number of people which can live in domains 
with given properties should be proportional to the total volume of these domains. If P p ((p, t\x) is 
stationary, then it seems reasonable to use it as a measure of the total volume of these domains at 
any given moment of time t. A detailed discussion of this distribution and all possible problems 



involved is contained in [|32|| . Some possible applications of this probability distribution for finding 
the most probable values of the effective gravitational constant in the Brans-Dicke inflationary 
cosmology can be found in [BvJ. We have found there that inflation in the Brans-Dicke theory 
leads to division of the Universe into different exponentially large domains with different values of 
the gravitational constant, and, correspondingly, with different values of density perturbations. 
We have shown also that the probability distribution P p {4>, t\x) in certain versions of this theory 
can be non-stationary, running away towards infinitely large 0, corresponding to infinitely large 
values of the effective Planck mass. To obtain a correct interpretation of these results one should 
develop some combination of stochastic approach to inflation and anthropic principle. 



A step in this direction was made in our recent paper with Juan Garcia-Bellido [Q, where 
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we investigated possible consequences of stochastic effects during inflation in application to the 
cosmological constant problem. We have shown, in particular, that by combining the baby 
universe theory with the Starobinsky model one can solve the problem of cosmological constant. 
The idea is that the overall growth of the total volume of all domains of the Universe containing 
matter with the density p ~ 10 -29 g-cm -3 is proportional to e Al ', where Ai is a parameter which 
depends on other constants of the theory, including the cosmological constant. As we have shown 
in [34], in the context of the Starobinsky model this rate increases as the cosmological constant 
becomes smaller. However, this solution works only if one believes into the baby universe theory 
(which allows one to compare universes with different values of the cosmological constant), and 
only if the cosmological constant cannot become negative. Thus, it cannot be considered as a 
final solution of the cosmological constant problem, but at least we have a completely new way 
to address this issue. 

In what follows I will briefly describe some nonperturbative effects which may lead to a 
considerable local deviation of density from the critical density of a flat universe. Results which 
I am going to discuss are very recent; they have been obtained in a collaboration with Dmitri 



Linde and Arthur Mezhlumian [35 |. 



Let us consider all parts of inflationary universe which contain a given field at a given 
moment of time t. One may wonder, what was the value of this field in those domains at the 
moment t — H^ 1 ? The answer is simple: One should add to the value of its classical drift 
A0 during the time H~ l , A0 = (pH^ 1 . One should also add the amplitude of a quantum jump 
5(f). The typical jump is given by 50 = =t|^- At the last stages of inflation this quantity is by 
many orders of magnitude smaller than A<f>. However, in which sense jumps ±^ are typical? 
If we consider any particular initial value of the field 0, then the typical jump from this point 
is indeed given by ±^-. However, if we are considering all domains with a given and trying 
to find all those domains from which the field could originate back in time, the answer may 
be quite different. Indeed, the total volume of all domains with a given field at any moment 

of time t strongly depends on 0: P p (4>) ~ A /~^~ Al ~ 1()8 , see eq. (f4T|). This means that the 
total volume of all domains which could jump towards the given field from the value + 50 

will be enhanced by a large additional factor ^ p p > ~^^ ~ (l + ^) V ^ 1 ■ On the other hand, the 
probability of large jumps 50 is suppressed by the Gaussian factor exp(— 2n ^2 ). One can easily 

verify that the product of these two factors has a sharp maximum at 50 = Ai0- and the width 
of this maximum is of the order J^. In other words, most of the domains of a given field are 
formed due to jumps which are greater than the "typical" ones by a factor Ai0 ± 0(1). 

Our part of the Universe in the inflationary scenario with V((p) = ^0 4 is formed at ~ 5 (in 
the units M p = 1), and the constant Ai ~ 2^/6ir ~ 8.68 for our choice of boundary conditions. 
This means that our part of the Universe should be created as a result of a jump down which 
is about Ai0 ~ 40 times greater than the standard jump. The standard jumps lead to density 
perturbations of the amplitude 5p ~ 5 • lCT 5 p c (in the normalization of ||). Thus, according 
to our nonperturbative analysis, we should live inside a region where density is smaller than 
the critical density by about 5p ~ 2 ■ lCT 3 p c . As we already mentioned, the probability of such 
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fluctuations should be suppressed by exp^— 2n J^ ) , which in our case gives the suppression factor 
~ exp(— 10 3 ). It is well known that exponentially suppressed perturbations typically give rise 
to spherically symmetric bubbles. Note also, that the Gaussian distribution suppressing the 
amplitude of the perturbations refers to the amplitude of a perturbation in its maximum. It is 
possible that we live not in the place corresponding to the maximum of the fluctuation. However, 
this could only happen if the nonperturbative jump down was even greater in the amplitude that 
we expected. Meanwhile, as we already mentioned, the distribution of the amplitudes of such 
jumps has width of only about J^. This means that we should live very close to the center of 
the giant fluctuation, and the difference of energy densities between the place where we live and 
the center of the "bubble" should be only about the same amplitude as the typical perturbative 
fluctuation Sp ~ 5 • 1CT 5 /9 C . In other words, we should live very close to the center of the nearly 
perfect spherically symmetric bubble, which contains matter with a smaller energy density than 
the matter outside it. 



It is very tempting to interpret this effect in such a way that the Universe around us becomes 
locally open, with 1 — Q ~ 10~ 3 . The true description of this effect is, of course, much more 
complicated, perhaps we should see the Hubble constant slightly growing at large distances. 
Note that in the simplest theory which we considered the magnitude of this effect is rather small. 
However, this smallness is proportional to the coefficient Ai, which is related to the maximal 
value of the Hubble constant during inflation. In some inflationary models based on the Brans- 



Dicke theories the effective value of Ai may be extremely large or even infinite [33, 34]. Thus, 



the typical deviation of density of the Universe from the average value in certain models can be 
very significant. If this is correct, then the degree of the local decrease of density can be used to 
determine the value of the parameter Ai. This means that by studying the large-scale structure 
of the Universe one can get information about the properties of space-time at Planckian densities! 



This effect is extremely unusual. We became partially satisfied by our understanding of 
this effect only after we confirmed its existence by four different methods, including computer 
simulations ||3~5| . However, there are still many problems associated with it. We are still trying to 



verify the existence of this effect and to find an adequate interpretation of our results. The most 
important problem is related to the choice of measure in inflationary cosmology. If we repeat our 
investigation using instead of P p (4>,t) the distribution P p (0, r), where r ~ log a, we do not get 
any nonperturbative effects of that kind. The source of this ambiguity can be easily understood. 
The total volume of all parts of an eternally inflating universe is infinite in the limit t — > oo (or 
t — > oo). Therefore when we are trying to compare volumes of domains with different properties, 
we are comparing infinities. This leads to answers depending on the way we are making this 



comparison, see the discussion of this issue in p2|-p 



Until we resolve the problem of measure in quantum cosmology, it will be hard to prove that the 
effect discussed above should actually take place. On the other hand, it is fair to say that now we 
cannot exclude the possibility of such an effect, and this by itself is a very unexpected conclusion. 
Few years ago we would say that such effects definitely contradict inflationary cosmology. Now 
we can only say that this is an open question to be studied experimentally, and that maybe 
cosmological observations can tell us something new about the nature of quantum gravity. 
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9 Conclusions 



The first realistic versions of inflationary cosmology have been proposed more than ten years ago. 
Since that time inflationary theory gradually became a very popular tool for investigation of the 
physical processes in the early universe. It also provides us with a good framework to study the 
large-scale structure of the Universe. However, inflationary theory is still at the stage of its active 
development. Within the last one or two years we have realized that inflation may end not only 
by slow rolling or by bubble formation, but in a kind of waterfall regime, associated with hybrid 
inflation. We have learned that reheating after inflation occurs in a much more complicated way 
than we thought before. I expect the theory of reheating to become a subject of intensive study 
in the near future. We have also learned that our understanding of the way inflation solves the 
monopole problem was absolutely incomplete, since we neglected the possibility that topological 
defects can exponentially grow. Finally, it was realized that inflationary theory may give rise to 
the theory of a stationary universe. This development completely changes our understanding of 
the structure and the fate of our Universe, and of our own place in the world. It is amazing that 
after many years of its active development inflationary theory still reveals so many surprising 
features and provides us with so many challenging problems to be solved! 
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